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Because  of  the  increasing  costs  of  resources,  the  increasing  attempts 
at  sabotage  or  destruction  of  these  resources  by  individuals  or  small  groups, 
whose  methods  have  become  more  sophisticated,  and  the  increasing  costs  of 
security  personnel,  the  development  of  intruder  detection  systems  has  become 
a  prime  objective  of  USAF,  of  DoD,  and  of  many  other  government  agencies.  One 
intruder  detection  concept,  which  has  many  possible  system  configurations,  uses 
a  leaky  coaxial  cable  as  the  primary  sensor.  To  properly  use  this  type  of 
cable,  the  properties  of  the  cable  should  be  understood  and  the  cable  performance 
known  under  a  wide  variety  of  conditions  and  situations.  The  work  in  this 
contract  deals  with  determining  analytically  and  experimentally  the  properties 
and  performance  of  leaky  coaxial  cables. 

Nuisance  alarms  plague  most,  if  not  all,  intruder  detection  systems.  A 
nuisance  alarm  is  one  generated  at  the  sensor  by  any  source  that  exhibits 
signatures  similar  to  those  of  personnel,  vehicles  or  craft  (e.g.  wildlife,  or 
fish).  The  ability  to  identify  an  intrusion  from  the  intruder  characteristic 
signature  would  decrease  the  nuisance  alarm  rate,  thereby  increasing  the 
ef f ect i veness  of  the  intrusion  system.  Results  indicate  that  it  is  possible, 
in  principle,  to  obtain  a  library  of  target  signatures,  and  from  these  determine 
criteria  for  identifying  an  intrusion. 

Leaky  coaxial  cables,  used  as  sensors,  lay  on  the  ground  or  are  buried.  In 
either  instance,  the  weather  affects  their  performance,  and  this,  in  turn, 
affects  vital  system  parameters  such  as  target  sensitivity,  cable  profile,  and 
false  alarm  rate.  Experiments  investigated  the  relationship  between  changes  in 


the  environmental  medium  and  changes  in  the  various  parameters  of  the  intrusion 
system.  Knowing  this  relationship  greatly  helps  in  designing  leaky  coaxial  cable 
sensor  systems  that  will  operate  effectively  in  all  kinds  of  weather. 

The  work  described- in  this  report  adds  significantly  to  the  knowledge 
that  is  needed  to  design  a  leaky  coaxial  cable  sensor  system  that  meets 
government  requirements. 

NICHOLAS  V.  KARAS 
Contract  Monitor 


Performance  Characteristics 
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Leaky  Coaxial  Cables 


1.1  Introduction 

This  is  the  final  report  required  under  USAF  contract 
No.  F19628-77-C-0249 .  It  deals,  in  the  main,  with  the  sub¬ 
ject  of  guided  radar  and  discusses  a  number  of  experiments, 
analytical  procedures  and  finally,  conclusions,  of  a  project 
designed  to  update  the  scientific  information  available  on 
this  topic. 

Queen's  University  has  been  in  the  forefront  of 
research  devoted  to  obstacle  detection  techniques  using 
guided  radar,  and  various  staff  members  have  played  an 
important  role  in  the  development  of  this  field.  A  biblio¬ 
graphy  containing  many  of  the  publications  of  interest  is 
provided  at  the  end  of  this  report. 

The  first  successful  attempt  to  detect  an  obstacle 
using  guided  radar  was  performed  by  Mackay  in  1973  and  a 
reliable  technique  to  perform  the  task  was  reported  by 
Penstone  and  Mackay  in  1974.  A  joint  effort  by  Queen's 
University  and  Computing  Devices  Company  resulted  in  the 
development  of  the  first  field  equipment  in  1977  under  a 
contract  with  the  Canadian  Penitentiary  Services.  Further 
development  of  the  equipment  together  with  extensive  field 
trials,  laboratory  experiments  and  analytical  studies  has 
continued  under  various  contracts  with  the  USAF,  the  Canadian 
National  Research  Council  and  other  government  agencies  in 
Canada  and  the  U.S.  This  report  adds  to  the  general  body 
of  knowledge  that  is  steadily  growing  in  this  area  and  pro¬ 
vides  new  insights  to  problems  of  target  detection  and 
identification. 

The  report  is  mainly  concerned  with  the  following 

aspects : 

(a)  Environmental  studies 

(b)  Velocity  measurements 

(c)  Target  identification,  classification  and  signature 

analysis 

(d)  Development  of  a  single-cable  technique  for 

relatively  short  perimeters. 

The  study  of  various  phenomena  related  to  guided 
radar  is  fraught  with  difficulties.  In  a  general  sense, 
a  reasonable  model  of  the  propagation  of  energy  along  the 
medium  (i.e.  via  leaky  coaxial  cable  and  through  the  immediate 
surroundings  of  the  cable)  has  been  hypothesized  by  Beal  [1,2] 
which  appears  to  match  observation  more  closely  than  other 
workers  [ 3 i .  Indeed,  the  results  of  well-controlled  labora- 
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tory  experiments  agree  very  well  with  theory  for  relatively 
short  lengths  of  cable  [  4  ].  There  has  been  some  difficulty, 
however,  in  trying  to  expand  the  theory  to  fully  meet  the 
practical  situation,  in  which  very  long  lines  are  used  and 
the  environment  is  relatively  uncontrollable.  The  problem, 
essentially,  is  not  that  the  theory  is  incorrect,  but  rather 
that  it  is  extremely  difficult  to  characterize  the  environ¬ 
ment  of  a  long  perimeter  through  various  soil  formations  and 
climatic  conditions.  General  statements  as  to  the  performance 
of  guided  radar  systems  in  various  environments  have  indeed 
been  made,  but  the  ability  to  reliably  verify  the  conclusions 
is  limited  by  the  means  available  to  accurately  identify 
the  environment  surrounding  a  buried  cable.  Continuous 
and  long-term  experimentation  is  therefore  necessary  to 
provide  much-needed  data  to  compare  with  hypothesis.  With 
this  in  mind,  therefore,  this  project  has  concentrated  in  the 
main  on  providing  a  good  data  base  for  further  study.  Many 
of  the  experiments  performed  required  considerable  sophisti¬ 
cation  of  equipment  and  deployment  of  personnel  and  were 
designed  to  provide  data  on  a  number  of  aspects  of  the  pro¬ 
ject;  results  for  each  particular  purpose  were  often  obtained 
by  processing  the  general  data  base  in  a  particular  manner. 

While  the  theory  of  propagation  via  guided  radar  has 
received  much  attention,  often  because  of  its  application  to 
the  communications  field,  little  is  yet  known  of  the  interaction 
of  a  target  with  the  detection  zone.  General  information,  con¬ 
cerning  the  reflected  power  expected  from  various  targets,  the 
influence  of  a  target  on  the  incident  field  and  the  region  of 
detection  around  the  cables  has  been  reported  in  the  past.f  5  ] 
However,  little  work  has  been  done  on  precisely  what  consti¬ 
tutes  a  target  for  guided  radar.  This  information  could  be 
very  useful  for  developing  target  classification  schemes  which 
would  allow  a  reduction  in  the  false  alarm  rate  in  a  deployed 
system.  A  start  has  been  made  on  this  problem  and  the  results 
are  presented  in  this  report. 

At  the  present  time,  a  number  of  field  models  of 
guided  radar  instruments  are  being  deployed  on  a  trial  basis. 

In  addition,  it  appears  likely  that  further  commercial  develop¬ 
ment  of  several  single-cable  techniques  will  continue  over 
the  next  few  years.  The  success  of  many  of  these  attempts 
to  develop  reliable  systems  for  various  purposes  will  still 
depend  on  a  proper  understanding  of  the  pehnomena  involved 
and  on  the  interaction  between  various  targets  and  the 
detection  procedure.  This  report  is  intended  to  assist 
these  attempts  in  providing  reliable  up-to-date  data  on 
several  important  items. 

1.2  Outline  of  Report 

Section  2  of  this  report  is  a  description  of  the 
guided  radar  research  facility  at  Queen's  University.  Included 
is  a  summary  of  the  characteristics  of  the  leaky  coaxial  cables 
used  in  the  experiments.  A  brief  description  is  given  of  the 
signal  processing  equipment  and  automated  data  logging  system. 


Much  more  detailed  information  on  the  system  is  available 
in  theses  by  Vinnins  [  6  ]  and  Patterson  [  7  ). 

Section  3  presents  the  results  of  the  experiments 
designed  to  assess  the  effect  of  changing  weather  conditions 
on  the  performance  of  a  guided  radar  system.  This  work  was 
carried  out  during  the  winter  and  early  spring  of  1979. 

The  work  presented  in  Section  4  demonstrates  that  a 
guided  radar  system  can  be  used  to  measure  the  velocity  of 
targets  within  the  detection  zone.  Experiments  are  also 
reported  in  which  the  velocities  of  two  targets  simultaneously 
within  the  detection  zone  were  independently  measured. 

Section  5  describes  target  classification  experi¬ 
ments  performed  using  a  variety  of  different  targets  including 
human  beings,  a  dog,  a  car  and  a  wooden  trolley.  An  explan¬ 
ation  is  given  of  some  of  the  distinctive  features  of  these 
signatures . 

One  approach  to  analysing  target  signatures  is  to 
simulate  the  expected  signal  return  from  targets  based  upon 
a  knowledge  of  their  motion  dynamics  and  the  characteristics 
of  the  guided  radar  system.  Some  preliminary  work  has  been 
done  using  this  approach  and  is  reported  in  Section  6. 

Section  7  deals  with  a  novel  single-cable  guided 
radar  suitable  for  relatively  short  perimeters.  A  technique 
is  described  which  allows  not  only  detection  but  also  location 
of  targets  even  though  the  system  operates  in  a  continuous 
wave  mode. 

Section  8  summarizes  the  results  obtained  during 
the  contract  period  and  also  contains  suggestions  for 
possible  future  research  in  guided  radar. 
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2.  System  Description 


2.1  Field  Site 

The  experimental  work  described  in  this  report  was 
carried  out  at  the  Queen's  University  test  site  situated 
approximately  eight  kilometers  northeast  of  Kingston.  The 
facility  consists  of  a  15  m  long  trailer  and  a  large  grass 
covered  field.  The  trailer  can  be  either  air  conditioned 
or  heated  to  maintain  the  inside  temperature  at  18°  C. 

A  plan  view  of  the  field  site  is  shown  in  figure  2.1. 
Along  the  western  edge  of  the  site  there  is  a  fence  line  and 
a  gravel  road.  Between  these,  the  ground  slopes  downward 
towards  the  road  and  is  tree  and  bush  covered.  Trees  also 
line  the  northern  boundary  of  the  site.  There  is  also  a 
65  metre  long  section  of  normal  gravel  railroad  bed  3  metres 
wide  and  30  cm  deep.  Standard  wooden  ties  and  tracks  are 
installed  on  the  gravel  bed  according  to  railroad  practice. 

The  leaky  coaxial  cables  were  normally  laid  on  the 
ground  or  suspended  at  various  heights  in  the  air  on  wooden 
stakes.  All  the  cable  insta l lations  were  in  the  open  part 
of  the  field  site  a  minimum  of  0  metres  away  from  any  trees. 

Figure  2.2  is  a  photograph  of  the  trailer  which 
housed  all  the  signal  processing  and  data  logging  equipment. 

A  section  of  the  railway  track  is  also  visible  in  the  photo¬ 
graph  as  is  the  wooden  moveable  platform  or  trolley  which 
was  used  as  a  target  in  some  of  the  experiments  described 
in  this  report.  The  trolley  was  propelled  by  a  rope  and 
pulley  arrangement  driven  by  a  3  h.p.  do  motor.  Figure  2.3 
is  a  photograph  of  the  central  cable  run  which  allowed  experi¬ 
ments  to  be  performed  using  a  variety  of  cable  spacings  and 
heights  above  the  ground.  Figure  2.4  is  a  photograph  taken 
inside  the  trailer  at  the  operator's  console  of  the  data 
logging  computer.  The  operator  had  a  good  view  of  the  test 
track  through  the  front  window  and  of  the  other  cable  runs 
through  the  side  window.  The  photograph  in  figure  2.5  shows 
the  analogue  equipment  which  made  up  the  guided  radar  system. 

2.2  Cables 

All  of  the  leaky  coaxial  cables  used  in  the  tests 
described  in  this  report  were  manufactured  by  Andrew  Antenna 
Company  Ltd.,  Ontario  Canaria  and  were  of  the  RADIAX  type. 

These  cables  are  constructed  in  an  identical  manner  to 
standard  HLLIAX  coaxial  cable  with  the  exception  that  slots 
have  been  milled  in  the  solid  outer  conductor  prior  to  the 
plastic  sheath  being  applied.  It  is  these  slots  that  cause 
the  cables  to  be  "leaky". 

A  family  of  RADIAX  cable  types  are  manufactured  by 
Andrews,  differing  only  in  the  size  of  the  slots  in  the 
outer  conductor.  These  have  been  named  by  the  manufacturer 
RX4-1,  RX4-2  and  RX4-3,  in  ascending  order  of  slot  size  and 
consequently  "leakiness".  All  three  cable  types  were  used  at 
the  field  site.  Manufacturer's  specifications  for  these  cables 
are  given  in  Table  2.1.  A  photograph  showing  the  cable 
construction  is  given  in  figure  2.6. 
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Plan  View  of  Field  Site 
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Table  2.1 

SPECIFICATIONS  FOR  LEAKY  CABLES  AT  60  MHz 


RX  4-1 

RX4-2 

RX4-3 

Nominal  size 

1/2" 

1/2" 

1/2" 

Impedance ,  ohms 

50 

50 

50 

Velocity,  percent  of 
light  in  free  space 

velocity  of 

79 

79 

79 

Typical  VSWR 

1.3 

1.3 

1.3 

Attenuation,  in  free 

space  (dB/100  ft) 

0.6 

0.75 

1.1 

Coupling  Loss  (dB  at 

20  ft.) 

83 

76 

71 

The  standard  cables  used  were  a  total  of  470  feet  in 
length  and  were  made  up  of  a  lead-in  section  200  ft.  long, 
350  ft.  of  radiating  cable  and  a  20  ft.  long  termination 
section.  The  lead-in  and  termination  sections  were  non¬ 
radiating  HELIAX .  The  centre  conductor  and  dielectric  were 
continuous  throughout  the  total  length  of  the  cables  and  no 
connectors  were  necessary  between  the  HELIAX  and  RADIAX 
portions.  N-type  connectors  were  installed  in  both  ends  of 
the  cables. 

The  cables  were  either  laid  on  the  surface  of  the 
ground  or  suspended  in  the  air  on  wooden  stakes.  A  variety 
of  cable  runs  were  installled  which  allowed  the  cables  to 
be  placed  at  different  heights  ranging  from  0  feet  to  6  feet 
and  at  different  spacings  from  5  feet  to  12  feet. 


2.3  Analog  Hardware 

2.3.1  Introduction 

The  guided  radar  transmitter  and  the  analog  portion 
of  the  receiver  is  built  up  of  commercially  available  pro¬ 
cessing  blocks  which  can  be  configured  as  either  a  pulsed 
or  continuous  wave  system.  The  interconnection  of  the 
various  units  is  made  through  coaxial  relays  which  can  be 
controlled  by  the  data  logging  computer.  This  allows  the 
system  to  operate  unattended  for  long  periods  of  time  with 
system  reconfiguration  occurring,  as  appropriate,  under  soft¬ 
ware  control. 

2.3.2  Pulse  Mode 

Short  bursts  of  an  RF  carrier  are  transmitted.  Three 
crystal  controlled  frequencies  are  available;  30,  60,  120  MHz. 
The  pulse  width  is  manually  controlled  and  was  normally  set 
to  60  ns  for  the  experiments  described  in  this  report.  The 
pulse  repetition  rate  is  controlled  by  the  PDP-11  computer 
and  was  normally  2.2  usee. 
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The  receiver  consists  of  a  broadband  preamplifier, 
followed  by  a  synchronous  detector  and  a  7  MHz  bandpass 
amplifier.  The  local  oscillator  could  be  phase  shifted 
by  90°  under  computer  control  allowing  either  the  in-phase  (I) 
or  quadrature  (Q)  component  of  the  return  signal  to  be 
monitored. 

A  block  diagram  of  the  pulse  system  is  shown  in 
figure  2.7. 


2.3.3  CW  Mode 

In  the  CW  mode  a  continuous  carrier  is  transmitted  at 
the  selected  frequency.  The  receiver  consists  of  a  hybrid, 
to  provide  the  required  quadrature  local  oscillator,  a 
pair  of  mixers,  a  pair  of  chopper  stablized  amplifiers  with 
switch  selected  gains  and  a  computer  controlled  electronic 
switch  to  select  which  received  signal  component  is  to  be 
monitored.  The  continuous  wave  system  is  shown  in  block  form 
in  figure  2.8. 

A  summary  of  some  of  the  important  characteristics 
of  the  analog  portion  of  the  system  is  given  in  Table  2.2. 


2.4  Digital  Hardware 

2.4.1  Introduction 

The  digital  portion  of  the  guided  radar  system  consists 
of  three  sections:  1)  Sampler  and  digital  preprocessor, 

2)  preset  and  timing  logic  and  3)  PDP-11/10  minicomputer. 

These  will  be  briefly  described  in  the  following  subsections. 

A  block  diagram  of  the  system  is  given  in  figure  2.9. 


2.4.2  Sampler  and  Digital  Preprocessor 

The  sampler  and  preprocessor  section  of  the  system 
performs  the  digitization  of  the  base-band  return  signal. 

Band-limited  components  of  the  return  signal  are 
multiplexed,  and  fed  by  the  R.F.  hardware  into  a  high  speed, 
sample  and  hold  module  which,  in  turn,  feeds  the  analog 
sampled  values  into  a  high  speed  8-bit  successive  approximation 
analog-to-digital  converter.  In  order  to  increase  the  signal- 
to-noise  ratio  of  the  return  signal,  many  successive  samples 
of  each  component  are  added  together.  Due  to  speed  limitations 
the  PDP  11/10  cannot  perform  this  operation  in  a  reasonable 
time  and  hence  a  preprocessor  unit  is  required  to  sum  con¬ 
secutive  samples  of  the  components  in  each  phase.  Upon  com¬ 
pletion  of  the  required  number  of  integrations,  the  integrated 
data  is  passed  on  to  the  computer,  ready  for  storage  on  the 
magnetic  tape.  Figure  2.10  is  a  block  diagram  of  the  organ¬ 
ization  of  the  preprocessor. 


Table  2.2 

GUIDED  RADAR  SYSTEM  SPECIFICATIONS 

General 

Carrier  Frequency 

Output  Power 
RF  Bandwidth 
Receiver  noise  figure 


40  ns  min,  manually  controlled 
450  KHz  max,  computer  controlled 
>  8  ft  in  RAD I AX 
7  MHz 

CW  Mode 


Pulse  Mode 

Pulse  width 

Pulse  repetition  rate 

Resolution 

Baseband  Bandwidth 


1  MHz  to  450  MHz  variable 
30,  60,  120  MHz  crystal 
controlled 

10  watts  peak 

50  MHz 

6  dB 


Baseband  Bandwidth 


32  Hz 


Figure  2.8  CW  Analog  System 


TO  PDP-II/IO 


Figure  2.10  Preprocessor 


2.4.3  Preset  and  Timing  Logic 

The  preset  and  timing  logic  section  controls  the 
events  within  the  system  in  accordance  with  the  system 
control  parameters.  The  timing  logic  works  on  an  internal 
50  MHz  TTL  clock  from  which  it  synchronizes  the  events  in  the 
system.  Thus,  on  command  from  the  timing  logic  the  selection 
of  either  the  in-phase  or  quadrature  phase  by  the  electronic 
switch,  the  issuing  of  the  sample  commands  to  the  sample  and 
hold,  and  the  start  conversion  to  the  ADC,  take  place,  in 
the  proper  sequence.  The  preset  logic  also  stores  the  required 
number  of  integrations  and  keeps  track  of  the  number  of  inte¬ 
grations  performed  by  the  preprocessor,  flagging  the  computer 
on  completion. 

2.4.4  The  PDP  11/10  Minicomputer 

The  minicomputer  serves  as  the  control  centre  and  per¬ 
forms  the  function  of  managing  the  operation  of  the  guided 
radar  system.  Operating  parameters  are  entered  through  a 
teletype  system.  Operating  parameters  are  entered  through  a 
teletype  terminal.  The  relevant  information,  in  turn,  is 
conveyed  to  the  different  sections  within  the  system.  Upon 
completion  of  all  operations,  including  the  preprocessing 
of  the  data,  the  computer  is  flagged  by  the  preprocessor  and 
reads  the  data  through  a  DR-11  digital  interface.  This 
data  can  be  processed  in  real  time  and  the  results  displayed 
on  an  oscilloscope  through  a  digital-to-analog  converter 
(DAC) ,  or  alternatively,  can  be  stored  on  DEC-tape  for  later 
analysis  on  a  larger  PDP- 15  computer  located  at  the  University. 


2.5  System  Parameters 

2.5.1  Introduction 

The  guided  radar  system  has  an  inherent  flexibility 
of  operation  wherein  the  adjustment  of  certain  system  para¬ 
meters  is  possible,  either  manually  or  by  program  control 
through  the  computer.  This  is  necessary  because  the  system 
is  designed  to  be  a  scientific  tool  rather  than  an  operational 
detection  system. 

The  relevant  parameters  of  the  system  are  the  following 

1.  Sampling  frequency  rate 

2.  Number  of  integrations 

3.  R.F.  carrier  frequency 

4.  Transmit  power. 


2.5.2  Sampling  Frequency  Rate 

The  sampling  frequency  rate  is  computer  controlled. 

The  doppler  frequency  in  the  return  signal  is  a  measure  of 
the  velocity  of  a  target.  Under  Nyquist's  criterion,  the 
sampling  frequency  rate  should  be  at  least  twice  the  maximum 
frequency  to  be  measured.  Hence  the  selection  of  the  sampling 
rate  will  be  determined  by  the  range  of  velocities  of  the 
targets  to  be  detected.  For  this  system,  the  rate  can  be 
varied  between  24  KHz  and  450  kHz.  As  will  be  seen  in  the 
next  section,  the  combination  of  this  parameter  and  the 
number  of  integrations  will  determine  the  fastest  detectable 
target  velocity. 

2.5.3  Number  of  Integrations 

To  increase  the  signal-to-noise  ratio  of  the  system, 
successive  sampled  values  are  integrated.  While  integration 
of  the  sampled  data  improves  the  signal-to-noise  ratio,  it 
also  decreases  the  effective  sampling  frequency  rate  of  the 
system. 

Thus,  in  any  experiment,  an  optimum  choice  of  the  two 
parameters  is  necessary.  The  number  of  integrations  performed 
on  each  component  can  be  selected  within  the  range  of  1  to 
4096. 

2.5.4  R.F.  Carrier  Frequency 

The  R.F.  carrier  frequency  can  be  manually  selected 
to  be  either  30  MHz,  60  MHz  or  120  MHz.  Several  parameters 
affect  the  particular  choice  of  operating  frequency.  The 
cables  have  a  useful  performance  range  of  20  MHz  to  400  MHz. 
Reduction  in  leakage  at  low  frequencies  and  attenuation 
effects  at  higher  frequencies  limit  their  usefulness  outside 
this  range.  In  the  study  of  moviny  targets,  the  choice  of  a 
particular  frequency  is  determined  by  the  maximum  expected 
target  velocity.  Since  the  doppler  frequency  shift  is  directly 
proportional  to  target  velocity  and  carrier  frequency,  the  use 
of  a  higher  carrier  frequency  provides  a  better  resolution 
in  the  detection  and  measurement  of  the  velocity  of  slower 
moving  targets,  while  a  lower  carrier  frequency  is  more  apt 
for  faster  targets  that  might  otherwise  have  been  outside 
the  system  bandwidth  as  determined  by  the  effective  sampling 
rate . 


2.5.5  Transmit  Power 

The  maximum  continuous  transmit  power  of  the  system 
is  5  watts.  However,  the  return  signal  power  varies  from 
experiment  to  experiment  depending  upon  factors  such  as  the 
size  of  the  slots  in  the  outer  conductor,  the  separation  between 
the  cables,  and  the  prevailing  environmental  and  weather  con¬ 
ditions.  Thus  for  each  experiment,  the  transmit  power  is 
chosen  to  be  the  highest  possible,  consistent  with  the  dynamic 
range  of  the  receiver. 


3  Environmental  Effects 


3.1  Introduction 

From  past  experience  it  is  known  that  the  type  of 
ground  upon  which  leaky  cables  are  laid  affects  their  per¬ 
formance  in  a  guided  radar  system.  Furthermore  the  condition 
of  the  ground  —  whether  it  is  frozen,  wet,  moist,  dry  —  also 
has  considerable  influence  on  system  parameters  such  as  target 
sensitivity,  cable  profile  and  false  alarm  rate.  The  purpose 
of  the  experiments  described  in  this  chapter  was  to  investi¬ 
gate  the  relationship  between  changes  in  the  environmental 
conditions  and  changes  in  various  parameters  of  a  guided 
radar  system.  The  work  was  carried  out  during  the  winter  and 
early  spring  seasons  when  the  most  significant  changes  in 
ground  conditions  would  be  expected  to  occur. 

An  automated  data  logging  system  was  set  up  which 
allowed  a  variety  of  measurements  to  be  made  on  both  a  pulse 
and  CW  guided  radar  system  at  regular  intervals.  The  details 
of  the  operation  of  the  system  are  given  in  Section  3.2  and 
the  measurements  made  are  described  in  Section  3.3.  The 
results  and  a  discussion  of  them  is  given  in  Sections  3.4  to  3.9 
Section  3 . 10  is  a  summary  of  this  phase  of  the  project. 


3.2  Data  Logging  Facility 

A  pair  of  RADIAX  RX4-2  cables,  350  feet  in  length 
were  used  for  the  monitoring  of  the  effect  of  weather  con¬ 
ditions  on  the  guided  radar  system.  These  cables  were  laid 
parallel  to  each  other,  approximately  6  feet  apart  along  the 
western  boundary  of  the  field  site.  Their  location  has  been 
marked  on  figure  2.1  ,  the  plan  view  of  the  site.  The 
cables  were  laid  on  the  surface  of  the  ground  several  months 
before  monitoring  began. 

Six  soil  moisture  cells  and  three  soil  temperature 
cells  were  also  installed  near  the  trailer  end  of  the 
cables.  A  detailed  schematic  of  the  cable  placement  is 
given  in  figure  3.1. 

The  computer  based  data  recording  facility  is  shown  in 
figure  3.2.  The  data  was  recorded  on  DECtape  on  the  PDP/11 
minicomputer  and  transferred  to  the  main  university  campus 
for  processing  on  a  PDP/15.  Mode  selection  hardware  consisting 
of  a  group  of  coaxial  relays  and  a  logic  block  allowed  either 
the  CW  hardware  or  the  pulse  system  hardware  to  be  selected 
under  software  control.  The  CW  and  pulse  hardware  were 
described  in  Section  2 . 3 . 3  and  2 . 3 . 2  respectively ,  of  this 
report.  A  computer  controlled  stepping  relay  allowed  the 
soil  moisture  and  temperature  cells  to  be  sequentially 
sampled  when  desired. 
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Figure  3.1:  Cable  Placement 


Figure  3.2:  Environmental  Data  Recording  Facility 
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3.3 


Measurements 


3.3.1  CW  Profile 

The  system  components  were  configured  as  shown  in 
figure  2.8  .  Each  time  a  CW  profile  measurement  was  made 
250  samples  of  both  the  I  and  Q  components  of  the  return 
signal  were  recorded.  Each  sample  of  a  particular  component 
was  the  sum  of  4096  individual  readings  and  the  time  taken 
for  measurement  was  4.5  seconds.  The  effective  profile 
sampling  rate  was  27  Hz. 


3.3.2  Pulse  Profile 

Figure  2.7  shows  the  pulse  guided  radar  system  used 
for  these  measurements.  A  window,  consisting  of  50  range 
cells,  a  total  of  107  metres  long  was  monitored.  The  window 
started  35  metres  from  the  transmitter  end  of  the  cables 
and  was  positioned  as  shown  in  figure  3.1. 

Five  complete  profiles  were  recorded  during  an  inter¬ 
val  of  4.5  seconds.  Therefore,  the  effective  sampling  rate 
was  1.1  Hz. 


3.3.3  Sensitivity 

The  system  sensitivity  was  measured  at  various  times 
by  means  of  a  standard  walk.  The  target,  a  human,  moved  down 
the  length  of  the  cables  midway  between  them  and  the  resulting 
target  response  was  used  as  an  indication  of  system  sensitivity. 
A  standard  walk  was  carried  out  every  three  or  four  days 
between  February  14th  and  April  1st.  More  frequent  walks 
were  taken  during  rapid  changes  in  environmental  conditions. 

3. 3. 3.1  Data  Format 

The  system  response  during  the  standard  walk  was 
recorded  using  the  pulsed  guided  radar  system  described  in 
Section  2.3.2. 

The  system  parameters  were: 

Pulse  width 
Repetition  rate 
Number  of  integrations 
Number  of  range  cells 
Range  cell  length 
Window  length 
Window  position 


60  ns 
456  kHz 
4096 
10 

8  ft. 

80  ft. 

350  ft.  from  transmit 
end  of  the  cables 


A  complete  profile  consisted  of  4096  x  10  samples 
of  each  return  signal  component  (I  and  Q) .  Therefore  the 
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effective  profile  sampling  rate  was  4096xlbx2  =  Hz* 

The  data  was  recorded  directly  onto  DECtape  at  the  field 
site  and  then  transferred  to  the  departmental  computer 
system  for  processing.  Figure  3.3  shows  the  format  in 
which  the  data  was  recorded  at  the  field  site. 


3. 3. 3. 2  Method  of  Walking 

The  return  signal  was  recorded  for  a  period  of 
approximately  34  seconds  before  the  target  entered  the 
detection  zone.  The  target  then  positioned  himself  between 
the  cables  at  the  point  marked  in  figure  3.1  and  proceeded 
to  walk  at  a  velocity  of  0.9  m/sec.  The  track  of  the  target 
was  approximately  midway  between  the  two  cables. 

The  walk  started  30  metres  before  the  observation 
window  and  continued  for  15  metres  past  the  window.  The  target 
stopped  at  marked  locations  every  27  m  for  a  period  of  five 
seconds . 

The  walks  were  done  under  a  variety  of  weather 
conditions;  consequently  there  were  variations  in  how  the 
target  moved  from  one  walk  to  another.  The  target  took 
small  steps  (<0.3  m)  if  possible,  to  minimize  the  influence 
of  leg  motion  on  the  results.  The  feet  remained  as  close 
to  the  ground  as  possible.  However,  when  there  were  large 
accumulations  of  snow  it  was  impossible  to  maintain  this 
style  of  walk.  Consequently,  a  number  of  different  walk 
styles  were  used. 

When  the  snow  was  deep  and  uncompacted  the  target 
still  took  small  steps  but  had  to  lift  his  feet  a  consider¬ 
able  distance  above  the  ground  for  each  step.  Since  all 
walks  were  taken  along  the  same  route  footprints  became  well 
defined  in  the  snow  and  were  used  as  guides  for  subsequent 
walks.  In  very  deep  powder  snow,  snowshoes  were  used. 

At  various  times  during  the  winter  a  thick  ice  crust 
formed  on  the  accumulated  snow  base.  This  crust  was  thick 
enough  to  support  the  weight  of  the  target  and,  therefore, 
the  target  walked  above  the  ground  at  a  height  equal  to  the 
snow  depth. 

Although  the  exact  conditions  of  the  target’s  walk 
varied  throughout  the  measurement  period,  the  differences 
in  style  would  not  be  expected  to  greatly  influence  the 
results. 


3.3.4  Environmental  Data 


The  soil  moisture  and  soil  temperature  cells  were 
sampled  each  time  an  automated  profile  measurement  was  made. 
The  resulting  data  was  recorded  on  DECtape. 

A  log  of  the  environmental  conditions  was  made  for 
each  visit  to  the  field  site.  The  data  recorded  included 
air  temperature,  general  weather  conditions  and  snow  depth 
at  six  selected  positions  along  the  length  of  the  cable. 

More  detailed  environmental  data  was  available  from 
the  Environment  Canada  weather  office  at  Norman  Rogers 
airport.  This  office  is  located  approximately  five  miles 
south-west  of  the  field  site.  Monthly  summaries  from 
Environment  Canada  contain  hourly  air  temperature  readings, 
average  snow  depth  on  the  ground,  wind  velocity  and  direction, 
information  concerning  the  onset  of  rain  or  snow  and  general 
weather  conditions 


3.4  General  weather  conditions 

The  majority  of  the  experiments  dealing  with  the 
effect  of  weather  conditions  on  the  leaky  cable  sensor  were 
performed  between  January  8,  1979  and  April  25,  1979.  This 
period  of  time  includes  the  portion  of  the  winter  season 
when  snowfall  is  heaviest  and  snow  cover  is  most  variable 
as  well  as  the  spring  thaw  and  early  spring  conditions. 

An  overview  of  the  weather  conditions  can  be  obtained  from 
the  records  maintained  at  Environment  Canada,  Kingston 
Weather  Office.  Based  upon  the  weather  summaries  kept 
there,  the  conditions  for  the  relevant  months  in  1979  were 
typical  of  those  occurring  in  other  years. 

Monthly  summaries,  for  January  to  April  1979,  pro¬ 
vided  by  Environment  Canada  for  the  Kingston  Airport  Weather 
Office,  are  attached  to  this  report  as  Appendix  A.  Relevant 
data  from  these  summaries  have  been  plotted  in  figures  3.4 
and  3.5.  Figure  3.4  is  a  plot  of  the  daily  maximum  and 
minimum  temperatures.  Also  plotted  is  the  daily  rainfall. 
Figure  3.5  is  a  plot  of  the  snow  on  the  ground  and  the  daily 
snowfall . 

The  temperature  ranged  from  -32°C  on  February  10  to 
+19°C  on  April  20.  There  were  several  periods  of  warm 
weather  during  the  winter  when  the  ground  was  snow  covered 
and  when  the  temperature  was  above  0°C  and  precipitation, 
in  the  form  of  rain,  occurred.  These  were  separated  by  cold 
spells.  Snow  on  the  ground  ranged  from  0  cm  to  45  cm 
with  the  largest  snowfalls  being  in  January.  The  snow  cover 
dropped  rapidly  at  the  beginning  of  March,  disappearing 
completely  on  March  13.  There  was  also  considerable  rain 
during  this  period.  The  end  of  March  and  beginning  of  April 
was  also  a  period  of  generally  warm  weather  and  heavy  rainfall. 


3.5 


Field  Site  Data 


Between  January  31  and  March  3  the  snow  depth  at 
six  positions  along  the  length  of  the  cables  was  measured 
each  time  a  standard  walk  was  performed.  These  results  are 
plotted  in  figure  3.6.  Due  to  drifting  of  the  snow  there  are 
variations  in  the  accumulated  snow  cover  from  one  position  to 
the  next.  The  average  snow  cover  at  the  field  site  is  com¬ 
pared  to  the  snow  depth  data  obtained  from  the  weather 
office  in  figure  3.7.  The  two  measurements  are  in  general 
agreement  indicating  that  weather  office  data  can  be  used 
when  average  field  site  data  is  unavailble  or  unreliable 
because  of  rapidly  changing  conditions. 

An  informal  log  was  kept  of  the  air  temperature  and 
the  general  weather  conditions  each  time  the  field  site  was 
visited. 

A  reading  was  made  of  the  soil  moisture  and  soil 
temperature  sensors  each  time  a  system  profile  measurement 
was  made.  The  sensors,  however,  had  been  installed  in  the 
ground  many  months  before  the  winter  freeze-up.  Unfor¬ 
tunately,  most  of  the  sensors  had  failed  during  the  early 
winter  and  could  not  be  removed  from  the  ground  to  be 
repaired  and  recalibrated.  Since  the  data  obtained  from 
these  sensors  could  not  be  considered  to  be  reliable  it 
will  not  be  referred  to  further  in  this  report. 


3.6  System  Noise 

The  effect  of  the  environmental  conditions  on  the 
system  noise  was  assessed  from  two  different  types  of 
data,  namely;  1)  the  initial  blocks  of  data  recorded  during 
standard  walks  when  no  target  was  within  the  sensor  detection 
zone,  and  2)  the  250  profile  samples  taken  during  each  CW 
profile  measurement. 

3.6.1  Standard  Walks 

During  the  recording  of  the  first  twelve  blocks  of  data 
for  each  standard  walk  no  target  was  present  within  the 
detection  zone  of  the  sensors.  These  data  consist  of  190 
samples  of  each  component  of  the  profile  for  each  of  the  10 
cells  recorded  over  a  period  of  41  seconds.  The  effective 
sampling  rate  was  1.1  Hz.  The  standard  deviation  of  these 
samples  is  an  indication  of  the  system  noise. 

Table  3.1  gives  the  standard  deviation  of  the 
magnitude  (M  )  of  the  profile  for  cells  1,  5  and  8  at  the 
time  of  eachbstandard  walk.  This  data  is  plotted  in  figure 
3.8  along  with  the  snow  depth  data.  The  noise  appears  to  be 
independent  of  the  snow  depth  and  the  weather  conditions. 

Figure  3.9  is  a  plot  of  Mgo  versus  cell  number  for 
three  different  dates.  A  typical  profile  is  also  plotted  in 
this  figure.  The  noise  from  a  particular  cell  is  not  related 
to  the  profile  magnitude  for  that  cell.  In  fact,  the  noise 
magnitude  appears  to  be  constant  throughout  the  length  of 
the  data  window. 
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Table  3.1 


SUMMARY  OF  NOISE  MEASUREMENTS 

DATE  TIME  MgD  (mV)  WEATHER  CONDITIONS 


Cell  1 

Cell  5 

Cell  8 

2  1 

5  p.m. 

3.9 

3.2 

5.0 

rain 

25 

12  p.m. 

4.4 

4.5 

4.4 

crusty  snow 

26 

4  |  > .  in . 

5.1 

4.8 

5.3 

heavy  snow  with  driftini! 

28 

5  p.m. 

4.0 

4.3 

4.1 

sunny 

3 

3 

p.m. 

3.6 

4.0 

3.8 

sunny 

5 

11 

a  ,m. 

3.9 

3.8 

4.8 

showers ,  snow  cover 
almost  gone 

6 

10 

a  .m. 

3.6 

3.8 

5.6 

rain 

7 

5 

p.m. 

3.8 

3.9 

4.9 

cloudy,  warm,  traces 
of  snow  on  ground 

10 

7 

a  .m. 

3.7 

4.9 

7.2 

raining 

12 

11 

a.m. 

4.3 

6.4 

7.4 

light  snow,  showers 

15 

5 

p.m. 

5.4 

5.7 

6.7 

clear 

2b 

4 

p.m. 

5.1 

6.0 

10 

snow,  rain  showers 

5 

4  p.m. 

5.1 

4.9 

5.6 

light  snow  overnight 

6 

5  p.m. 

4.2 

4.6 

7.9 

heavy  snow  overnight. 

stormy  conditions. 

light  cover  of  snow 

9 

5  p.m. 

4.3 

5.0 

4.8 

snowing  all  day 

10 

1  a.m. 

3.9 

5.3 

6.1 

snow 

10 

3  p.m. 

3.7 

5.2 

6.4 

13 

3  p.m. 

4. 1 

3.5 

5.7 

overcast 

15 

2  p.m. 

18 

17 

21 

cloudy 

17 

4  p.m. 

5.7 

4.6 

5.8 

sunny 

Figure  3.6  Snow  Depth  Measured  at  Field  Site 


Figure  3.7  Snow  Depth  Comparison 
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Figure  3.8  MgD  versus  Date 
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Figure  3.9  M  versus  Cell  NuMber 
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Typical  noise  plots  are  given  in  figures  3.10 
and  3.11.  Figure  3.10  is  a  plot  of  I,  Q,  and  M  for  cells 
1  and  8  on  February  25.  Figure  3.11  is  a  plot  of  M  for 
cells  1  and  8  on  three  different  dates. 


3.6.2  CW  Profile 

Each  time  a  CW  profile  was  recorded  250  readings 
of  I  and  Q  were  made  at  a  profile  sampling  rate  of  27  Hz. 
The  total  time  taken  for  each  record  was  4.5  seconds.  The 
standard  deviation  of  the  250  samples  was  calculated. 

From  this  data  the  noise  on  the  CW  profile  appeared 
to  be  unrelated  to  the  weather  conditions,  snow  depth  or 
ground  conditions. 


3.6.3  Receiver  Noise 

Two  types  of  receiver  noise  tests  were  performed. 

In  the  first  type  the  noise  generated  within  the  receiver 
itself  was  measured  by  terminating  the  input  with  a  standard 
50Q  load  and  recording  the  CW  profile  once  per  hour  for  66 
hours.  In  the  second  test  the  cable  was  left  connected  to 
the  receiver  input  but  the  transmitter  was  turned  off.  The 
CW  profile  was  again  measured  once  per  hour  for  66  hours. 

In  this  case  the  data  includes  not  only  receiver  noise  but 
also  any  noise  in  the  environment  picked  up  by  the  receive 
sensor. 

The  data  was  processed  in  the  standard  fashion  to 

yield 

I  -  the  average  of  250  samples  of  the  in-phase  component 

Q  -  the  average  of  250  samples  of  the  quadrature  component 

M  -  the  average  magnitude 

ISD  “  the  standard  deviation  of  the  250  I  samples 

Qsd  ~  the  standard  deviation  of  the  250  Q  samples 

Figure  3.12  is  a  plot  of  I,  Q  and  M  with  the  receiver 
input  terminated  and  figure  3.13  is  a  plot  of  the  results  with 
the  receive  cable  connected.  They  are  essentially  the  same. 

The  stability  of  the  Q  channel  is  considerably  poorer  than 
that  of  the  I  channel. 

The  maximum  value  of  Igp  and  Qsd  was  the  same  whether 
or  not  the  receive  cable  was  connected  indicating  that  the 
majority  of  the  noise  was  generated  within  the  receiver 
rather  than  being  caused  by  external  sources  of  electro¬ 
magnetic  energy. 

A  summary  of  the  receiver  noise  tests  is  given  in 
Table  3.2. 
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Table  3.2 

RECEIVER  NOISE  TESTS 

Receive  Cable  Receiver  Input 

Connected  Terminated 


STABILITY  (mV) 


j  .  t> .  *4  Conclusions 

The  majority  of  the  noise  within  the  system  is 
generated  in  the  receiver  itself.  The  environmental  con¬ 
ditions  do  not  add  appreciably  to  the  high  frequency  noise. 
At  the  Oueen’s  University  field  site  externally  generated 
electromagnetic  fields  do  not  significantly  increase  the 
noise  level.  Since  the  noise  is  receiver  generated,  it  is 
uniform  in  magnitude  along  the  length  of  the  sensors  and  is 
not  related  to  the  profile  or  sensitivity  at  a  particular 
location . 


3.7  Sensitivity 

3.7.1  Typical  Results 

Table  3.3  is  a  summary  of  the  date  and  time  when 
each  standard  walk  was  performed.  Included  in  the  table 
is  a  brief  description  of  the  weather  conditions  at  the  time. 

Figure  3.14  is  a  plot  of  the  I  component  of  the 
profile  during  a  typical  standard  walk  for  each  of  the  ten 
range  cells;  figure  3.15  is  a  plot  of  the  corresponding  mag¬ 
nitude  response  for  each  cell  during  the  same  walk.  In  these 
plots  the  abscissa  variable  is  time  and  the  ordinate  variable 
is  target  response.  Each  curve  is  the  response  from  a 
different  range  cell.  The  profile  magnitude  is  also  plotted 
on  this  figure  for  reference.  During  this  walk  there  was 
11  cm  of  snow  on  the  ground.  Although  the  snow  had  a  thick 
ice  crust  on  top  the  footprints  from  previous  walks  were 
distinct.  Therefore,  the  target  lifted  his  feet  a  con¬ 
siderable  distance  above  the  ground  during  each  step.  This 
is  the  probable  cause  of  the  rapid  variations  seen  in  the 
sensitivity  plots. 

An  alternate  presentation  of  the  data  recorded  during 
the  same  standard  walk  is  given  in  figure  3.16.  In  this 
figure  target  response  is  plotted  versus  the  cell  number 
(distance  along  the  cable)  at  equal  time  intervals  during 
the  walk.  The  width  of  the  target  response  remains  con¬ 
stant  as  the  target  walks  through  the  data  "window" 
implying  that  the  distance  resolution  of  the  system  also 
remains  constant.  This  fact  was  true  for  all  standard  walks, 
independent  of  the  weather  conditions. 


3.7.2  Sensitivity  Variation  with  Weather  Conditions 

The  maximum  magnitude  of  the  target  response  in 
each  cell  was  used  as  a  measure  of  system  sensitivity.  The 
sensitivity  for  each  cell  during  each  walk,  relative  to  the 
sensitivity  for  the  corresponding  cell  on  the  reference  date 
of  February  25,  is  plotted  in  figure  3.17.  The  curves  for 
all  ten  cells  are  plotted  on  the  same  ordinate  scale.  How¬ 
ever  for  clarity  they  have  been  shifted  vertically  relative 
to  each  other.  The  absolute  sensitivity  of  each  cell  is 
marked  at  the  beginning  of  each  curve.  The  rainfall  data 
from  the  weather  office  is  also  plotted  in  this  figure. 
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Table  3.3 


SUMMARY  OF  STANDARD  WALKS 


DATE 

TIME 

SNOW  DEPTH 
(cm) 

TEMPERATURE 

(cm) 

WEATHER  CONDITIONS 

— 

Feb.  21 

5  p.m. 

35 

+  1 

rain 

25 

12  p.m. 

11 

crusty  snow 

26 

4  p.m. 

11 

heavy  snow  with  drifting 

28 

5  p.m. 

18 

+  3 

sunny 

— * 

Mar.  3 

3  p.m. 

6 

+  6 

sunny 

5 

11  a.m. 

2 

+  3 

showers ,  snow  cover 

almost  gone 

6 

10  a.m. 

2 

+  3 

rain 

7 

5  p.m. 

2 

+  3 

cloudy,  warm,  traces 

of  snow  on  ground 

10 

7  a.m. 

2 

-1 

raining 

12 

11  a.m. 

1 

0 

light  snow,  showers 

15 

5  p.m. 

0 

-10 

clear 

26 

4  p.m. 

0 

-3 

snow,  rain  showers 

Apr.  5 

4  p.m. 

2 

+  5 

light  snow  overnight 

6 

5  p.m. 

1 

0 

heavy  snow  overnight, 

stormy  conditions, 

light  cover  of  snow 

9 

5  p.m. 

0 

0 

snowing  all  day 

10 

1  a.m. 

5 

-1 

snow 

10 

3  p.m. 

5 

+  16 

13 

3  p.m. 

0 

+  13 

overcast 

— 

15 

2  p.m. 

0 

+  10 

cloudy 

17 

4  p.m. 

0 

+  15 

sunny 
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Figure  3.16  System  Response  -  Standard  Walk 
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Figure  3.17  Sensitivity  Variation  -  Each  Cell 


The  sensitivity  varies  in  a  similar  fashion  for  each 
range  cell  in  response  to  changes  in  the  weather  conditions. 
Moisture  either  on  the  ground  or  saturating  the  snow  above 
the  cables  has  the  greatest  effect  on  sensitivity. 

Between  February  21  and  February  24  a  total  of  20.8  mm 
of  rain  fell  and  the  snow  depth  decreased  from  35  cm  to  15  cm. 
By  the  25th  the  snow  was  well  saturated  with  moisture.  The 
sensitivity  decreased  by  5  to  7  dB  during  this  period  of  time. 
Between  March  2  and  March  6  rain  and  warm  temperatures  caused 
the  rapid  disappearance  of  the  remaining  snow  cover.  The 
sensitivity  was  reduced  by  a  further  5  dB.  From  the  7th  of 
March  to  the  14th  of  March  the  ground  was  clear  of  snow  and 
was  gradually  drying  out  although  the  surface  was  still  moist. 
The  sensitivity  rose  by  approximately  10  dB  to  a  maximum  on 
March  15th.  More  complete  drying  of  the  surface  after  March 
15th  resulted  in  a  gradual  decline  in  sensitivity. 


3.7.3  Sensitivity  Variation  from  Cell  to  Cell 

The  variation  in  sensitivity  in  response  to  changing 
weather  conditions  along  the  length  of  the  cable  was  relatively 
uniform.  This  is  demonstrated  by  figure  3.18.  For  this  plot, 
the  sensitivity  on  February  25th  was  used  as  a  reference  for 
each  cell.  Only  the  change  in  sensitivity  from  the  February 
25th  reference  value  for  each  cell  is  plotted  in  figure  3.18. 
The  greatest  difference  in  the  change  in  sensitivity  occurred 
on  April  17  between  cell  1  and  cell  10.  The  difference  was 
6.2  dB. 

Table  3.4  is  a  summary  of  the  sensitivity  changes 
during  the  experimental  period.  The  change  in  sensitivity  at 
a  particular  point  on  the  cable  appears  to  be  unrelated  to 
the  absolute  sensitivity  at  that  point. 


3.7.4  Summary 

The  sensitivity  is  most  clearly  influenced  by 
moisture  in  the  vicinity  of  the  cables,  either  in  the 
surface  soil  beneath  the  cables  or  in  the  form  of  moist  snow 
above  the  cables.  Under  very  wet  conditions  (rapidly  melting 
snow)  a  large  decrease  in  sensitivity  was  noted  averaging 
9.1  dB .  With  the  ground  moist  but  clear  of  snow  and  surface 
water  (March  16)  the  sensitivity  appeared  to  be  2  or  3  dB 
greater  than  when  the  ground  was  dry.  Changes  in  sensitivity 
were  relatively  uniform  along  the  length  of  the  observation 
window . 

3.8  Pulse  Profile 

Two  types  of  measurement  were  made  of  pulse  profile. 

The  first  series  of  measurements  was  made  immediately  prior 
to  each  standard  walk.  These  results  will  be  discussed  in 
sub  sections  3.8.1  to  3.8.3.  Pulse  profile  was  also  recorded 
at  hourly  intervals  for  certain  periods  of  time  (See  section 
3.3.2).  These  results,  along  with  the  corresponding  CW  profile 
will  be  discussed  in  section  3.9. 


Table  3.4 


VARIATION  IN  SENSITIVY 


RANGE 

CELL 

MAX  (dB) 

SENSITIVITY 

MIN  (dB)  VARIATION 

ABSOLUTE 
SENSITIVITY 
(dB)  Feb.  25 

1 

9.5 

-4.1 

13.6 

0 

2 

8.9 

-4.5 

13.4 

.63 

3 

4.1 

-3.5 

7.6 

5.4 

4 

4.2 

-2.6 

6.8 

11.8 

5 

5.0 

-3.4 

8.4 

14.0 

6 

5.2 

-3.9 

9.1 

16.7 

7 

5.2 

-4.1 

9.3 

18.3 

8 

7.8 

-4.9 

12.7 

16.4 

9 

7.6 

-4.3 

11.9 

15.2 

10 

9.3 

-6.5 

15.8 

16.7 

Note: 

All  sensitivities 

are  calculated  relative 

to  the  sensitivity 

for  a  particular 

range  cell  on 

February  25 

. 
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3.8.1  Profile  Variation  with  Weather  Conditions 

The  recordings  of  pulse  profile  were  made  immediately 
prior  to  each  standard  walk.  Figure  3.19  is  a  plot  of  the 
profile  in  each  cell  with  the  profile  on  February  21st  being 
chosen  as  the  reference  for  each  cell.  The  rainfall  data 
is  also  plotted  in  this  figure.  The  three  symbols  x  .  and  ° 
are  used  in  the  plots  of  the  response  of  the  different  range 
cells  merely  for  clarity. 

The  pattern  of  profile  variation  is  similar  for  all 
cells  during  most  of  the  measurement  period.  The  largest 
profile  changes  occurred  in  response  to  rain  and  melting 
snow. 


3.8.2  Profile  Variation  from  Cell  to  Cell 

Figure  3.20  is  a  plot  of  the  change  in  profile  along 
the  cable  on  various  dates  when  standard  walks  were  performed 
The  profile  in  each  cell  has  been  normalized  to  the  corres¬ 
ponding  profile  on  February  21st. 

The  change  in  profile  is  relatively  uniform  along 
the  cable.  Cells  #4  and  #8  have  anomolously  low  profiles 
on  some  of  the  dates  when  an  experiment  was  performed.  The 
reason  for  this  is  unknown  although  it  could  result  from 
very  localized  differences  in  ground  moisture  (puddle  for¬ 
mation  perhaps,  caused  by  a  local  ground  depression) . 

Table  3.5  summarizes  the  changes  in  profile  in 
each  range  cell.  Also  given  is  the  reference  profile  on 
February  21,  normalized  to  cell  1. 


3.8.3  Sensitivity  versus  cell  profile 

A  question  which  has  arisen  in  the  past  concerning 
guided  radar  is  whether  sensitivity  and  profile  are  directly 
related  in  a  particular  range  cell.  A  comparison  of 
figures  3.17  and  3.19  indicates  that  both  these  parameters 
respond  in  a  similar  fashion  to  changes  in  the  environmental 
conditions  although  there  does  not  appear  to  be  a  direct 
correlation  between  them.  This  is  further  illustrated  by 
figure  3.21  in  which  profile  is  plotted  versus  sensitivity 
for  cell  1  and  cell  5. 


3.9  Automated  Profile  Monitoring 

3.9.1  Results 

At  various  times  between  December  21,  1978  and 
April  9,  1979  the  guided  radar  system  was  configured  for  auto¬ 
mated  sampling  of  CW  profile  and  pulse  profile.  Measurements 
of  the  profile  were  made  at  a  pre-selected  rate,  usually  once 
per  hour.  At  times  more  frequent  measurements  were  made. 

Table  3.6  is  a  summary  of  the  dates  and  weather  conditions 
during  these  measurements. 


Table  3 . 5 


PULSE  PROFILE  VARIATION 


PROFILE  MAGNITUDE 


RANGE 

CELL 

MAX  <dB) 

MIN  <dB) 

VARIATION  (dB) 

1 

7.7 

-  .92 

8.6 

2 

6.1 

-  .45 

6.6 

3 

6.5 

-1.5 

8.0 

4 

8.2 

-14 

22.2 

5 

6.4 

-  6.8 

13.2 

6 

6.3 

-  2.9 

9.2 

7 

8.1 

-  .29 

8.4 

8 

5.2 

-16 

21.2 

9 

3.7 

-11 

14.7 

10 

5.5 

-  7.5 

13.0 

REFERENCE 
PROFILE 
Feb.  21 

0 

2.5 
4.8 
6.3 

11.3 

13.4 

9.5 
11.3 
16.0 
15.1 


Note:  All  profiles  are  calculated  relative  to  the  profile  for  a 
particular  range  cell  on  February  21. 


Figure  3.21:  Sensitivity  versus  Profile 


Table  3.6 


SUMMARY  OF  AUTOMATED  PROFILE  RECORDINGS 


START 


Dec . 

21 

12:49 

p.m 

Jan. 

6 

1:05 

a  .m 

10 

3:  30 

p.m 

17 

3:55 

a  .m 

19 

6:15 

p.m 

20 

9:15 

a.m 

27 

6:28 

p.m 

30 

5:42 

a.m 

30 

3:25 

p.m 

30 

9:38 

p.m 

Feb. 

1 

7:34 

a.m 

3 

1:02 

p.m 

6 

6:31 

a.m 

9 

12:51 

p.m 

12 

1:15 

a.m 

16 

9:30 

p.m 

23 

8:00 

p.m 

28 

6:30 

a.m 

Mar . 

7 

9:30 

a.m 

8 

10:40 

p.m 

14 

10:26 

p.m 

18 

4:36 

p.m 

Apr . 

2 

1:05 

a  .m 

6 

1:00 

a.m 

9 

5:19 

p.m 

LENGTH  RATE 
(hours)  (per  hour) 


66 

1 

51 

1 

26 

1 

45 

1 

8 

1 

33 

1 

33 

1 

20 

10 

33 

10 

20 

1 

33 

1 

66 

1 

33 

1 

55 

1 

64 

1 

66 

1 

33 

1 

50 

1 

50 

■> 

*• 

33 

1 

33 

1 

33 

2 

19 

1 

66 

1 

7 

2 

WEATHER 


sunny,  cold 
overcast,  cold 
very  cold,  heavy  snow 
very  cold,  overcast 
rain 

warm,  overcast 
cold,  overcast 
light  snow 
light  snow 

clear,  sunny,  cold 
clear,  cold 
clear,  very  cold 
clear,  very  cold 
overcast,  cold 
light  dry  snow 
warm,  heavy  rain 
clear,  warm 

overcast,  warm 
overcast,  warm 
rain 

clear,  warm 

warm,  heavy  rain 
rain  and  snow 
heavy  wet  snow 
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The  large  amount  of  data  can  most  conveniently 
be  displayed  graphically.  Figure  3.22  is  an  example  of  a 
plot  of  the  pulse  profile  each  hour  for  a  period  of  33 
hours  starting  at  9  p.m.  on  December  22.  What  is  of 
interest  in  this  data  is  the  change  in  profile  in  response 
to  changes  in  the  weather.  To  make  the  changes  more  evident 
the  first  profile  in  the  sequence  can  be  selected  as  a 
reference  and  only  the  difference  between  a  profile  and 
the  reference  plotted.  This  has  been  done  to  produce 
figure  3.23.  The  weather  conditions  are  marked  on  the 
plots.  Figure  3.24  is  the  CW  profile  (I,  0  and  M)  for  the 
same  period  of  time. 

At  the  beginning  of  this  time  period  the  ground  was 
frozen  at  the  surface  and  there  was  a  trace  of  snow  on  the 
ground.  Overnight  temperatures  reached  a  low  of  -4°C. 

The  profile  was  stable  until  the  morning  of  the  23rd.  Tem¬ 
peratures  of  +2°C  caused  the  snow  to  disappear  and  the  ground 
surface  to  melt.  The  CW  profile  magnitude  dropped  by  2.1  dB 
due  to  the  ground  moisture.  The  profile  rose  back  to  its 
former  level  when  the  temperature  again  fell  below  the 
freezing  point.  This  effect  was  seen  several  times  during 
December . 

Figures  3.25  and  3.26  show  the  pulse  and  CW  profile 
changes  during  clear  stable  weather.  The  temperature  was 
below  -10°C  throughout  and  19  cm  of  snow  covered  the  cables. 
The  maximum  CW  profile  variation  was  0.7  dB. 

Rain  on  previously  dry  ground  had  a  large  effect  on 
both  the  CW  and  pulse  profiles.  This  is  illustrated  in 
figures  3.27  and  3.28.  Approximately  1  mm  of  rain  and  5  cm 
of  snow  fell  overnight  on  the  6th  of  April  with  the  sky 
clearing  in  the  early  morning.  The  CW  profile  rose  by  9  dB 
during  the  precipitation  and  dropped  by  6.2  dB  as  the  ground 
dried  during  the  warm  sunny  day.  Rain  during  the  evening 
of  the  7th  caused  the  profile  to  rise  again. 

The  effect  of  rain  saturating  previously  dry  snow 
is  shown  in  figures  3.29  and  3.30.  Initially , in  these 
figures, a  light  drizzle  was  falling.  This  changed  to 
moderate  rain  around  noon  of  the  23rd  followed  by  gradual 
clearing.  Approximately  20  cm  of  snow  was  on  the  ground  on 
the  23rd  which  decreased  to  15  cm  during  the  rainfall. 

There  is  a  delay  between  the  onset  of  the  rain  and  an  effect 
being  seen  in  the  profiles  as  might  be  expected.  Once  the 
moisture  had  percolated  through  the  snow  cover  to  the 
vicinity  of  the  cables  significant  changes  were  evident  in 
both  the  pulse  and  CW  profiles. 

3.10  Summary 

The  weather  conditions  obviously  significantly 
affect  the  performance  of  a  guided  radar  system.  From  the 
data  collected  it  is  possible  to  draw  general  conclusions 
as  to  the  nature  of  these  effects.  However,  since  the 
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Figure  3.24:  CW  Profile  -  December  22 
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Figure  3.30:  CW  Profile  -  rain  on  snow 


monitoring  was  done  throughout  only  one  winter  season 
and  since  a  particular  combination  of  weather  conditions 
occurs  only  once  or  twice  during  a  season,  in  an  uncontrolled 
fashion,  it  is  not  possible  to  be  precise  about  the  relation 
between  the  magnitude  of  an  effect.  The  following  general 
conclusions  can  be  made. 

The  presence  of  large  amounts  of  moisture  either  in 
the  snow  above  the  cables  or  on  the  ground  beneath  has  the 
most  significant  effect  on  profile  and  sensitivity.  A  decrease 
of  up  to  10  dB  in  the  sensitivity  can  occur  in  any  particular 
range  cell.  Similar  variations  in  the  profile  in  any  range 
cell  can  also  be  expected. 

The  sensitivity  is  enhanced  by  2  to  3  dB  when  the 
ground  is  only  slightly  moist.  The  pulse  profile  also 
increases  under  the  same  conditions. 

The  change  in  sensitivity  in  any  range  cell  is 
unrelated  to  the  absolute  sensitivity  in  the  same  cell. 

Profile  variations  of  up  to  22  dB  were  measured 
during  the  experimental  period  although  10  dB  was  a  more 
normal  figure. 

The  profile  and  sensitivity  in  a  particular  range 
cell  appear  to  be  only  loosely  correlated. 

Rapid  changes  in  CW  profile  can  be  expected  when 
rain  falls  on  previously  dry  ground.  An  increase  of  9  dB 
at  a  rate  of  0.6  dB/hr .  was  noted.  A  decrease  of  1  dB/hr 
was  seen  when  the  ground  was  drying  out. 

Thawing  of  the  ground  surface  causes  a  decrease  of 
about  2  dB  in  the  CW  profile. 

Noise  in  the  system  is  generated  principally  within 
the  receiver .  Weather  conditions  do  not  affect  the  noise 
appreciably. 


4 .  VELOCITY  MEASUREMENT 


4.1  Introduction 

A  target  moving  longitudinally  between  the  two  sensors 
in  a  ported  coaxial  cable  intrusion  detection  system  causes 
a  doppler  shift  in  the  frequency  of  the  return  signal.  The 
purpose  of  the  experiments  described  in  this  section  was  to 
demonstrate  that  the  magnitude  of  the  doppler  shift  could 
be  used  to  measure  the  target  velocity  and  furthermore,  that 
the  direction  of  motion  of  the  target  with  respect  to  the 
transmitter  could  also  be  determined  from  the  return  signal. 

The  principle  upon  which  the  measurement  technique  is 
based  can  be  conveniently  described  with  the  aid  of  a  vector 
model  of  the  return  signal  in  the  presence  of  a  target  shown 
in  figure  4.1.  Synchronous  detection  is  used  to  obtain  the 
in-phase  (I)  and  quadrature  (Q)  components  of  the  return 
signal.  The  reference  profile,  the  return  signal  in  the 
absence  of  a  target,  is  represented  in  the  figure  by  vector 
R.  With  a  target  present,  the_profile  changes  to  vector  P  . 
The  vector  difference  between  R  and  P  is_ref erred  to  as 
the  target  profile  T.  The  phase  <J>(t)of  T  is  calculated 
from  the  relation 


<Mt)  =  tan  1  ( AQ/  AI ) 

where 


and 


AQ  = 


°R 


(4.1) 


The  magnitude  of  T  is 


1/2 

M-  =  [  (AI)  2  +  (AQ)  2  ] 


(4.2) 


As  a  target  moves  down  the  cable  the  target  profile  T 
rotates  about  the  centre  of  the  circle  in  figure  4.1  with  an 
angular  velocity  which  is  directly  proportional  to  the  linear 
velocity  of  the  target. 

The  magnitude  of  the  target  profile  T  is  determined  by 
the  magnitude  of  the  field  in  the  detection  zone  and  the 
cross-section  of  the  target.  The  cross-section  is  assumed 
to  be  constant  for  a  particular  target  in  a  particular 
orientation.  The  phase  of  T  depends  upon  the  time  taken 
for  the  electromagnetic  energy  to  propagate  from  the  trans¬ 
mitter  to  the  target  and  back  again.  As  the  target  moves 
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CABLE  PROFILE  -  TARGET  PRESENT 


IN  PHASE  COMPONENT 


Figure  4.1 


Vector  Model  of  the  Return  Signal 


away  from  the  transmit  end  of  the  cables  this  path  length 
increases  and  hence,  so  does  the  phase  of  T.  As  the 
target  moves  towards  the  transmitter  the  path  length  and 
consequently,  the  phase  of  T  decreases. 

A  360  degree_phase  shift,  or  one  complete  rotation 
of  target  profile  T  in  the  vector  diagram,  corresponds  to 
the  target  moving  a  distance  of  half  the  wavelength  of  the 
RF  signal  within  the  leaky  co-axial  cable.  If  a  target 
moves  at  a  constant  velocity, _T  rotates  at  a  constant 
angular  velocity,  and  AI  and  AQ  (the  projections  of  T 
on  the  orthogonal  axes)  will  vary  sinusoidally.  The  frequency 
of  these  sinusoids  will  be  directly  proportional  to  the  velo¬ 
city  of  the  target,  and  the  direction  of  rotation  of  T  will 
be  an  indication  of  the  direction  of  motion  of  the  target. 


4.2  Processing 

The  overall  processing  scheme  for  the  measurement  of 
velocity  from  the  return  signal  is  shown  in  Figure  4.2.  The 
recorded  experimental  data  consists  of  two  discrete  sampled 
time  sequences  corresponding  to  the  in-phase  and  quadrature 
components  of  the  return  signal.  Let  these  two  sequences 
be  {ln}  and  { Q  > ,  where 

1  =  sample  in  the  in-phase  sampled  time  sequence 

1  of  the  return  signal 

0n  =  nth  sample  in  the  quadrature  sampled  time  sequence 
of  the  return  signal 

n  =  10,  1,  2,  .  .  .,  N- 1 1 

N  =  number  of  samples  taken  in  each  array. 

The  reference  cable  profile  is  first  removed  from  the 
two  arrays .  If: 

IR  =  in-phase  component  of  the  reference  cable  profile, 
and, 


The  result  of  the  transform  operation  is  a  complex 
spectrum  F{x(n)),  where 


N-l 


N-l 


=  l  ( AI  +  j  AQn)exp(-j  k)  (4.6) 

_ A  *  * 


If  f  is  the  sampling  frequency,  k  takes  on  discrete 
values  that  are  integer  multiples  of  fs/N.  The  value  of 
f S/N  is  the  frequency  resolution  of  the  spectrum.  The 
magnitudes  of  the  elements  of  the  complex  array  form  the 
power  spectrum  P(k),  where 


N-l  ,_n 

P  (k)  =  £  (AI  +  j  AQn)exp(-j  k)  (4.7) 

n=0  n  w 


Using  this  power  spectrum,  the  average  velocity  v  of 
the  target  can  be  calculated  from  the  frequency  of  the  largest 
spectral  component  by  the  relation: 

V  =  J  X  ff  x  (4.8) 

2  t  f  „ 


where. 


=  frequency  of  the  principal  spectral  component  in 
the  spectrum 

v  =  vc’ocity  of  propagation  of  the  electromagnetic 

>,oax  j  tj10  leaky  cable 

f  =  carrier  frequency 


4.3  Experimental  Procedure 

A  series  of  experiments  was  performed  to  verify  the 
analysis  presented  in  the  previous  section  of  this  report. 
In  each  experiment  the  target  moved  along  the  cables  at 
a  velocity  that  was  roughly  constant.  The  total  distance 
covered  by  the  target  and  the  time  taken  were  each  measured 
and  the  average  velocity  was  calculated. 
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The  return  signal  from  a  variety  of  different  targets, 
moving  at  different  speeds  in  both  directions  was  recorded. 

The  target  moved  either  midway  between  the  cables  or  just 
outside  the  receive  cable.  In  all  the  experiments  the  sampling 
rate  was  high  enough  and  the  target  velocity  was  low  enough 
that  no  aliasing  occurrred. 

The  cables  were  either  suspended  in  the  air  at  approxi¬ 
mately  1.5  metres  or  laid  on  the  ground.  Targets  used 
included  a  human  walking,  running  or  bicycling,  a  dog  trotting, 
a  car  and  a  wooden  cart.  The  cart  travelled  on  a  standard 
railroad  track  and  bed  which  had  been  installed  at  the  Queen’s 
University  field  site  for  other  experiments.  It  was  propelled 
by  a  D.C.  motor  through  a  rope  and  pulley  system. 

The  normal  sequence  of  events  in  each  experiment  was 
as  follows: 

1)  the  return  signal  was  recorded  for  5  seconds  with  no 
target  in  the  detection  zone.  This  was  done  to  establish 
a  reference  profile 

2)  the  target  entered  the  detection  zone  and  remained  at  the 
starting  point  for  5  seconds 

3)  the  target  moved  down  the  cables  at  the  desired  rate 

Obviously,  this  procedure  was  impractical  when  either 
the  trolley  or  the  dog  were  used  as  targets.  In  the 
experiments  involving  either  of  these  targets  an  estimate  of 
the  reference  profile  was  obtained  by  averaging  all  the  samples 
in  each  component  array  during  the  whole  experiment. 


4.4  Return  Signal 

Figure  4.3  is  a  plot  of  the  I  and  Q  components  of  the 
target  profile  recorded  while  a  target  (the  trolley)  moved 
between  the  cables  at  a  velocity  of  6.3  metres/sec.  For 
this  experiment  the  cables  were  in  the  air  on  wooden  poles 
just  outside  the  railway  tracks.  This  signal  is  typical  of 
those  recorded  in  all  the  experiments  dealing  with  velocity 
measurement . 

An  inspection  of  each  component  of  the  profile 
reveals  the  following  features: 


a)  a  basic  sinusoidal  siynal 

b)  a  low  frequency  overall  envelope 

c)  a  gradual  decrease  in  the  peak  value  of  the  envelope 

d)  a  small  amount  of  high  frequency  noise. 

The  frequency  of  the  basic  sinusoid  is  proportional 
to  the  velocity  of  the  target.  The  slowly  changing  magni¬ 
tude  envelope  results  from  non-uniform  sensitivity  of  the 
system  at  different  locations  along  the  path  due  to  inter¬ 
ference  between  two  different  modes  of  propagation  of  energy 
along  the  sensor  f 16  1 .  It  transforms  to  a  low  frequency 
in  the  spectrum  and  does  not  hinder  the  velocity  measurement 
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process.  The  gradual  decrease  in  the  peak  value  of  the 
envelope  as  the  target  moves  away  from  the  transmit  end  of  the 
sensor  is  due  to  cable  attenuation.  The  high  frequency  noise 
can  be  attributed  to  the  cables  vibrating  in  the  wind,  grass 
in  the  vicinity  of  the  cables,  noise  in  the  receiver  and 
quantization  noise. 


4.5  Velocity  Measurement 

Table  4.1  summarizes  the  different  combinations  of 
target,  target  speed  and  cable  placement  used  in  these 
experiments . 

Figure  4.4  shows  the  magnitude  of  the  power  spectra 
obtained  in  two  different  experiments.  In  each  case  the 
trolley  moved  along  the  track  at  a  speed  of  0.3  metres/sec. 
but  in  different  directions.  For  the  purposes  of  this  report, 
all  velocities  toward  the  transmit/receive  end  of  the  cables 
will  be  considered  to  be  positive,  while  all  velocities  away 
from  the  transmit/receive  end  are  negative.  Since  the 
sampling  rate  of  the  system  is  52.8  Hz,  the  power  spectrum 
covers  a  range  of  +  26.4  Hz.  However,  only  the  portion  of 
the  spectrum  that  contains  the  major  frequency  components  of 
the  baseband  signal  is  plotted.  A  total  of  1024  samples  of 
the  return  signal  were  used  in  the  transform  implying  that 
the  sample  points  in  the  transform  are  .052  Hz  apart. 

Since,  in  all  experiments,  the  velocities  of  the  targets 
were  fairly  constant,  the  longitudinal  axis  can  be  considered 
to  be  a  linear  measure  of  doppler  frequency  or  velocity,  with 
a  range  of  ±  0.5  Hz  or  f  1.0  metres/sec.  The  vertical  axis 
is  the  magnitude  of  the  complex  fourier  transform.  Although 
dimensionless,  this  magnitude  is  directly  proportional  to 
the  energy  at  each  discrete  frequency  in  the  baseband  signal. 
In  Figure  4.4  the  frequencies  of  the  largest  peaks  in  the 
power  spectra  correspond  to  the  doppler  frequencies  of  the 
moving  targets.  The  velocities  are  computed  using  equation 
4.8. 

The  velocities  computed  from  the  spectra  of  Figure  4.4 
were  found  to  be  in  agreement  with  the  measured  values, 
indicating  that  the  technique  was  effective  for  targets 
moving  in  both  directions. 

The  power  spectra  in  Figure  4.5  are  computed  from  the 
return  signals  of  five  experiments  in  which  a  human  walked 
down  the  cable  at  different  speeds  between  0.5  metres/sec. 
and  2.5  metres/sec.  at  intervals  of  0.5  metre/sec. 

The  graph  of  Figure  4.6  shows  the  relationship  between 
the  velocities  computed  from  the  spectra  of  Figure  4.5  and 
the  corresponding  velocities  measured  during  the  experiment. 

It  is  seen  that  over  this  range  of  velocities,  the  computed 
and  measured  values  are  in  good  agreement. 

The  spectrum  shown  in  Figure  4.7  was  obtained  from  a 
human  running  at  an  average  velocity  of  6.6  m/sec.  The  cables 
were  laid  on  the  ground.  The  main  spectral  component  occurs 
at  a  frequency  indicating  the  average  velocity  of  the  target. 
As  well,  secondary  peaks  occur  on  either  side  of  the  main 
component.  These  secondary  peaks  are  believed  to  be  related 
to  the  motion  of  the  legs  and  will  be  discussed  in  detail 
in  the  next  chapter.  73 
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Figure  4.7:  Spectrum  for  a  Human  Running 


Figure  4.8:  Spectrum  for  a  Human  Bicycling 


Figure  4.8  is  the  spectrum  for  a  human  bicycling,  and 
Figure  4.9  for  a  dog  running  between  the  cables.  Figure  4.10 
is  the  result  obtained  for  a  car  travelling  down  the  sensor 
just  outside  the  cables. 

In  all  cases  the  target  velocity  calculated  from  equation 
4.8  corresponded  to  the  measured  target  velocity.  Neither 
the  position  of  the  cables  with  respect  to  the  ground  nor  the 
location  of  the  target  with  respect  to  the  cables  influenced 
the  accuracy  of  the  velocity  measurement. 


4.6  Multiple  Targets 

A  number  of  experiments  were  conducted  in  which  two 
targets  moved  along  the  sensor  simultaneously. 

To  obtain  the  results  shown  in  Figure  4.11  a  human 
walked  at  +  0.85  m/s ec  in  front  of  the  trolley  which  was 
moving  in  the  same  direction  at  0.2  m/sec.  Spectral  com¬ 
ponents  are  clearly  visible  corresponding  to  the  average 
velocity  of  each  target.  In  a  different  experiment  two 
humans  walked  at  the  same  speed  in  opposite  directions.  The 
spectrum  of  the  return  signal  for  the  situation  is  shown  in 
Figure  4.12.  Again,  a  peak  occurs  in  the  spectrum  corresponding 
to  each  target. 


4.7  Conclusions 

A  C.W.  guided  radar  system  using  leaky  coaxial  cables 
as  the  sensor  can  be  used  to  measure  the  velocity  of  targets 
within  the  detection  zone  by  calculating  the  doppler  shift 
in  the  return  signal.  In  this  analysis  the  doppler  frequency 
was  calculated  by  means  of  a  discrete  complex  fourier  transform 
of  a  sequence  of  samples  of  the  return  signal.  If,  however, 
only  a  single  target  is  in  the  detection  zone  and  only  a  measure¬ 
ment  of  velocity  is  required  (rather  than  "target  signature" 
discussed  in  the  next  chapter)  a  much  simpler  processing 
scheme  is  possible.  Either  the  I  or  Q  components  of  the 
target  profile  T  could  be  passed  through  an  infinite  clipper 
followed  by  a  zero  crossing  detector.  The  zero  crossing 
frequency  would  be  directly  related  to  target  velocity. 

The  velocities  of  multiple  targets  within  the  detection 
zone  can  be  independently  measured  using  a  C.W.  system  if 
they  differ  from  each  other  (either  in  speed  or  direction) 
by  an  amount  greater  than  the  resolution  of  the  system  in 
frequency  space.  This  property  would  allow  the  detection  of 
multiple  targets  with  a  C.W.  system,  assuming,  of  course, 
that  their  velocities  differ.  If  the  targets  are  separated 
in  space  but  have  the  same  velocity  then  a  pulsed  guided 
radar  must  be  used  to  resolve  the  targets.  However,  if  the 
targets  are  resolvable  in  space  with  the  pulsed  system  then 
the  velocities  of  each  target  can  be  independently  measured. 
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Figure  4.9:  Spectrum  for  a  Dog  Running 
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Figure  4.10:  Spectrum  for  a  Car 


5. 


Target  Classif ication 


5.1  Introduction 

After  transforming  the  return  signal  from  a  leaky  coaxial 
sensor  into  frequency  space  it  is  possible  to  identify  com¬ 
ponents  in  the  calculated  spectrum  which  result  from  the 
doppler  shift  due  to  the  average  target  velocity.  The 
demonstration  of  this  fact  has  been  described  in  the  previous 
chapter.  The  spectrum,  however,  should  contain  more  infor¬ 
mation  than  simply  the  average  velocity  of  a  target;  it 
should  also  contain  information  about  the  characteristic 
motion  of  individual  parts  of  the  target  if  they  move  at  a 
velocity  other  than  the  average  velocity.  For  instance,  in  the 
case  of  a  human  walking,  although  the  body  moves  at  a  relatively 
constant  velocity,  the  arms  and  legs  do  not.  In  addition,  the 
body  and  each  of  the  appendages  have  different  effective  radar 
cross-sections  that  are  determined  by  their  individual  shape, 
size  and  location  relative  to  the  sensors.  The  overall  return 
signal  will  be  some  combination  of  the  return  signal  from  each 
body  part  and  hence  will  reflect  the  motion  characteristics 
of  the  target  components.  The  spectrum  should  also  contain 
detailed  information  about  target  characteristics . 

The  purpose  of  the  experiments  described  in  this  chapter 
was  to  demonstrate  that  the  spectrum  of  the  return  signal  con¬ 
tains  information  which  could  be  useful  for  target  classification. 


5.2  Processing 

The  processing  required  for  the  analysis  of  the  return 
signal  for  this  application  is  an  extension  of  that  used  in  the 
measurement  of  target  velocity.  The  scheme  is  shown  in  block 
form  in  figure  5.1.  As  was  described  in  section  4.1,  the 
recorded  data  consists  of  the  arrays  iln^  and  iOnK  (the  in- 
phase  and  quadrature  components  of  the  return  signal) ,  from 
which  the  target  vector  arrays  (A  I,.)  and  iAOn)  are  generated. 

For  performing  a  discrete  fourier  transform  the  finite 
duration  sampled  time  sequence  is  assumed  to  be  periodic  with 
a  period  equal  to  the  observation  interval.  The  transformation 
then  yields  a  discrete  spectrum.  However,  when  the  data  recorded 
consists  of  a  non-integer  number  of  cycles  of  a  quasi-per iodic 
waveform,  phase  jumps  in-  discontinuities  are  caused  at  either 
end  of  the  time  sequence  when  the  time  sequence  is  assumed  to 
be  periodic.  Spurious  spectral  component s  are  produced  by 
these  discontinuities.  The  application  of  a  digital  "window" 
to  the  time  sequences  before  transformation  aids  in  ridding 
the  frequency  domain  spectrum  of  some  of  this  spectral  leakage 
and  smooths  the  spectral  samples  in  the  fourier  transform. 

Extensive  literature  exists  concerning  the  selection 
of  an  appropriate  filter  or  window  for  particular  applications 
:  8,9].  Harris  ;  10  1  provides  a  good  summary  of  the  use  of 

various  windows  tor  harmonic  analysis.  Two  important  criteria 
for  the  selection  of  a  window  are  the  detectability  of  single 
t requency  tones  and  the  resolvability  of  neighbouring  tones  of 
similar  strength.  The  6  dB  bandwidth  of  the  central  lobe  of 
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Figure  5.1:  Target  Signature  Processing  Scheme 
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the  transform  of  the  window  is  a  measure  of  the  resolvability. 

A  low  sidelobe  level  increases  the  detectability  of  weak  tones 
in  the  presence  of  strong  ones.  Choice  of  a  window  usually 
entails  a  trade-off  between  the  highest  sidelobe  level  and 
the  6  dB  bandwidth. 

The  window  employed  in  the  processing  scheme  for 
target  classification  was  the  so-called  Blackman-Harris  window. 
The  weighting  functions  wn  for  this  window  are  generated  by 
the  relation: 


a0#  a 2#  a^  are  constants  selected  in  accordance 

with  desired  window  properties. 

The  Blackman  Harris  window  was  selected  because  it 
possessed  properties  suitable  for  use  in  harmonic  analysis, 
and  in  addition  was  simple  to  implement  digitally.  In 
the  trade-off  between  the  sidelobe  level  and  the  6  dB 
bandwidth,  a  sidelobe  level  of  -  74  dB  was  chosen,  commen¬ 
surate  with  the  dynamic  range  of  the  system.  The  coefficients 
of  the  window  were: 

aQ  =  0.40217 

ax  =  0.49703 

a2  =  0.09392 

a3  =  0.00183 


The  window  had  the  following  properties: 
highest  sidelobe  level  =  -  74  dB 
sidelobe  fall  off  rate  =  -  6  dB/octave 
6.0  dB  bandwidth  -  2.44  bins1 
3.0  dB  bandwidth  =  1.74  bins. 

in  a  discrete  spectrum  the  distance  between  adjacent  discrete 
frequencies  is  called  a  bin.  ete 


Figure  5.2(a)  shows  the  envelope  of  the  window,  and 
Figure  5.2(b)  the  log  magnitude  of  its  transform  in  dB, 
normalized  to  its  peak  value. 

The  data  arrays  (AIn)  and  {AQ^}  are  windowed  using 
this  Blackman-Harris  window  and  modified  to  arrays  {wn  AIn) 
and  (wn  AQn}  respectively. 

The  windowed  time  sequences  are  the  real  and  imaginary 
inputs  to  a  discrete  complex  fourier  transform.  If  the  complex 
input  sequence  is  (x(n)},  then: 


•jN-1 

4In  +  5  «„  }n=0 


(5.2) 


These  sequences  are  then  transformed  into  complex  frequency 
space.  Hence: 


F 


I_  +  j  w 
n  n 


*n 


} 


N-l 

n=0 


'  K  tun  4In  +  3  wn  Mln>“*,(-jTr  k) 


(5.3) 


where  N  is  the  total  number  of  samples. 

fs  is  the  sampling  frequency,  and  k  takes  on  discrete  frequency 
values  that  are  integer  multiples  of  fs/N. 

The  power  spectrum  P(k)  is: 


P(k) 


l  lwn  41n  +  3  wn  4°„>exp(-3Tr  k> 


(5.4) 


The  base  10  logarithm  of  this  power  spectrum 
generates  the  log  spectrum  referred  to  as  the  target  signa¬ 
ture.  Hence  the  target  signature  T(k)  is: 


T (k)  =  20  Iog1Q 


J  „(wn  AIn  +  3  wn  «Vexp(-3TT  k) 


(5.5) 
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Figure  5.2:  Blackman-Harris  Window 


This  target  signature  is  the  principal, tool  that  is 
used  to  demonstrate  the  feasibility  of  rudimentary  target 
classification  using  a  guided  radar  system. 


5.3  Experiments 

A  variety  of  different  types  of  targets  were  used  in 
the  experiments.  Since  the  purpose  was  to  demonstrate  that 
certain  characteristics  of  the  way  in  which  different  targets 
moved  would  be  reflected  in  the  target  signature,  targets  were 
chosen  with  kinematics  of  locomotion  as  different  from  each 
other  as  possible.  These  included  humans  walking,  running 
and  bicycling;  a  dog  running;  the  trolley  on  a  railroad 
track;  and  a  car  being  driven  along  the  sensor  just  outside 
the  cables. 

A  series  of  experiments  was  also  performed  using  a 
human  target  with  the  cables  at  various  heights  above  the 
ground  from  0  m  to  2  m.  These  experiments  were  designed  to 
demonstrate  the  variation  in  target  signature  which  occurs 
because  of  the  non  uniform  distribution  of  the  field  around 
the  leaky  cable  sensor. 


5.4  Results 

5.4.1  Human  Targets,  normal  walk 

The  first  series  of  experiments  all  involved  a  human 
as  the  target,  walking  midway  between  the  cables.  The 
walk  however,  was  not  a  natural  one.  To  reduce  the  complexity 
of  the  observed  signature  the  target  maintained  his  arms 
stationary  with  respect  to  his  body  at  all  times. 

The  experimental  series  is  summarized  in  Table  5.1  and 
typical  signatures  are  given  in  figures  5.3  to  5.5. 

Figure  5.3  demonstrates  the  effect  of  cable  height 
on  the  sidelobe  structure  of  the  signatures.  The  walk  which 
resulted  in  figure  5.3(a)  was  done  at  -  1  m/sec  while  the 
other  two  walks  were  at  -  0.5  m/sec.  As  might  be  expected 
a  large  spectral  component  occurs  in  all  signatures  corres¬ 
ponding  to  the  average  velocity  of  the  target.  There  is  also 
considerable  structure  in  the  signatures  at  other  frequencies 
particularly  when  the  cables  were  on  the  ground.  In  parti¬ 
cular,  in  figure  5.3(a)  it  is  seen  that  the  primary  peak  is 
located  at  -0.5  Hz,  while  secondary  peaks  occur  at  -4.0  Hz, 
-2.2  Hz,  +0.5  Hz,  +1.1  Hz,  +2.9  Hz  and  +4.6  Hz.  For  the 
present,  if  the  peak  at  +0.5  Hz  is  ignored,  then  all  other 
peaks  are  found  to  occur  at  frequencies  that  are  integer 
multiples  of  the  primary  frequency  component.  In  addition, 
the  distance  (in  frequency  space)  between  neighbouring  peaks 
is  approximately  three  times  the  primary  frequency.  The  first 

sidelobe  was  typically  13  dB  lower  than  the  central  peak  and 
the  second  side  lobes  were  approximately  38  dB  lower. 
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Table  5.1 


SUMMARY  OF  HUMAN  TARGET  SIGNATURE  EXPERIMENTS 


Speed  (m/s) 


Direction  Cable  Height  (m) 


0.5 

1.0 

1.5 
2.0 

2.5 
0.5 

1.5 
0.5 


0 

0 


0 

0 

0 

0 

0 

0 

1 

1 

2 

2 

0  Marking  Time 
0  Swinging  arms 
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(  VELOCITY  m/s  ) 


(  VELOCITY  m/s  ) 


The  signature  for  the  human  walk  with  the  cables  1  m 
in  the  air  is  shown  in  figure  5.3(b).  Again  sidelobes  are 
evident  in  the  signatures  although  they  are  considerably 
less  distinct  than  they  were  with  the  cables  on  the  ground. 

When  the  cables  were  2  m  in  the  air  (figure  5.3(c))  the 
signature  was  basically  similar  to  that  obtained  with  the 
cables  1  m  high.  In  both  these  figures  the  position  of  the 
discernable  sidelobes  is  related  to  the  position  of  the  main 
peak  in  the  same  fashion  as  was  noted  above. 

The  existence  of  the  secondary  peaks  in  the  signa¬ 
tures  for  a  human  walk  and  their  particular  position  in  fre¬ 
quency  space  with  respect  to  the  main  peak  are  due  to  the 
motion  of  the  legs  and  feet  of  the  target  and  the  cyclic 
manner  with  which  they  move. 

When  the  cables  are  on  the  ground  the  effects  of  leg 
motion  are  more  clearly  seen  in  the  signature  because  the 
coupled  field  strength  in  the  vicinity  of  the  legs  increases. 

The  energy  seen  in  the  return  signal  at  any  given  doppler 
frequency  would  be  expected  to  be  the  product  of  an  effective 
"cross-section"  for  the  portion  of  the  target  concerned  and 
the  field  strength  in  that  region. 

The  position  of  the  sidelobes  with  respect  to  the 
central  peak  is  determined  by  the  normal  pattern  of  leg 
motion.  For  a  walking  human  target,  these  secondary  peaks 
appear  in  the  vicinity  of  a  frequency 

f  =  f  +  3nf  5.6 

no  o 

where  fQ  =  frequency  representing  the  average  walking  speed 
n  =  an  integer 

This  relation  holds  true  over  a  considerable  range  of  walking 
speeds.  This  is  demonstrated  by  the  signatures  shown  in 
figure  5.4  and  the  data  from  a  series  of  experiments  summarized 
in  Table  5.2.  The  validity  of  equation  5.6  reflects  the 
fact  that  increased  walking  speed  is  normally  accomplished 
by  humans  by  increasing  step  frequency  rather  than  step 
length. 

The  return  signal  and  consequent  signature  for  a 
human  target  walking  between  the  sensors  has  been  simulated 
based  upon  a  study  of  the  kinematics  of  leg  motion  during 
normal  locomotion.  The  relationship  between  sidelobe  position 
in  the  signature  and  the  peak  representing  average  velocity  will 
be  discussed  more  fully  in  the  Section  6  of  this  report. 

In  addition  to  the  secondary  peaks  in  the  signature 
which  have  been  identified  as  resulting  from  leg  motion,  for 
each  principal  spectral  component  an  image  component  occurs, 

15  to  25  dB  lower,  at  the  corresponding  frequency  on  the 
opposite  side  of  the  zero  frequency  axis.  This  effect  is 
clearly  visible  in  the  signatures  shown  in  figure  5.5.  These 
peaks  occur  because  of  inaccuracies  in  the  quadrature  oscillator 
and  in  the  I  and  Q  channel  gains  in  the  synchronous  detector. 
They  will  be  more  fully  discussed  in  Section  5.5. 


5.4.2  Human  Target,  other  motions 

Signatures  were  also  obtained  for  human  targets 

a)  marking  time  (figure  5.6)  and  b)  swinging  his  arms  parallel 
to  the  cables  (figure  5.7).  In  the  first  case  the  cables  were 
on  the  ground  and  in  the  second  case  were  1  metre  off  the 
ground.  In  both  cases  the  cyclic  motion  of  the  appropriate 
appendage  produces  a  distinct  sideband  structure.  Since  the 
average  velocity  is  0  m/sec  the  largest  spectral  component 
occurs  at  zero  Hz. 


5.4.3  Other  T  arget s 

Four  other  targets  with  motion  dynamics  considerably 
different  from  those  of  a  human  walking  were  also  used.  These 
were 

a)  a  car  -  figure  5.8 

b)  the  trolley  on  the  railroad  track  -  figure  5.9 

c)  a  human  bicycling  -  figure  5.10 

d)  a  dog  running  between  the  cables  -  figure  5.11 

For  both  the  car  and  the  trolley  the  target  signatures 
were  composed  of  the  main  spectral  peak  corresponding  to  the 
average  velocity  of  the  target  with  few  other  distinct  features 
in  the  sidebands.  The  "leakage"  of  energy  from  one  frequency 
to  the  corresponding  frequency  on  the  opposite  side  of  the 
zero  frequency  axis  is  clear  in  these  signatures. 

The  signatures  for  the  bicycle  target,  shown  in 
figure  5.10,  have  some  structure  in  the  sidebands  but  the 
regularity  of  the  sidelobes  that  occurred  in  the  signatures 
of  walking  targets  is  not  present.  Although  there  is  some 
motion  of  the  legs  in  the  horizontal  plane  when  riding  a 
bicycle,  the  speed  of  travel  of  the  legs,  with  respect  to  the 
average  velocity,  is  less  than  when  walking  due  to  the  gearing  on 
the  bicycle.  Also,  the  amplitude  of  motion  in  the  horizontal 
plane,  again  with  respect  to  the  bicycle,  is  also  less  then 
when  walking. 

Figure  5.11  shows  the  target  signature  of  a  dog  running 
down  the  cables.  Here,  as  in  the  signature  of  the  human  walk, 
secondary  peaks  are  visible,  and  are  caused  by  the  motion  of 
the  dog's  feet.  In  frequency  space,  the  position  of  the 
secondary  peaks  relative  to  the  principal  peak  differs  from 
that  of  the  human  target  signature.  This  is  due  to  the  fact 
that  dogs  have  different  kinematics  of  motion  than  humans. 


5.5  System  Inaccuracies 

Certain  inaccuracies  in  the  system  hardware  components 
cause  the  spectral  components  in  the  signatures  to  be  reflected 
about  the  zero  frequency  axis.  Measurements  were  performed  on  the 
system  components  to  establish  the  magnitude  of  various  system 
errors  on  the  target  signature.  These  errors  were: 
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1.  Deviation  from  quadrature;  The  deviation  from 
quadrature  of  the  ANZAC  JH-10-4  broadband,  90  degree  hybrid 
was  measured  with  a  vector  voltmeter.  It  was  found  to  be  3.0 
degrees,  which  is  within  the  worst  case  specifications  at 

60  MHz,  as  quoted  by  the  manufacturer. 

2.  Amplitude  imbalance:  The  difference  in  gain  between 
the  I  channel  and  the  Q  channel,  measured  at  60  MHz  was  1  dB. 


The  influence  of  these  errors  on  the  target  signature 
can  be  demonstrated  analytically  for  the  simple  case  of  the 
return  signal  consisting  of  a  single  spectral  component. 
Errors  in  amplitude  or  in  phase,  change  a  spectral  component 
by  transferring  some  of  the  energy  at  that  frequency  to  its 
image  frequency. 

It  will  be  assumed  for  the  purposes  of  this  analysis 
that  the  observation  interval  is  of  duration  equal  to  the 
time  required  by  the  target  to  move  a  distance  corresponding 
to  half  the  wavelength  of  the  RF  signal  in  the  coaxial  cable. 
This  corresponds  to  a  baseband  return  signal  of  a  single 
cycle  in  each  component.  The  fourier  transform  of  this 
signal  places  the  corresponding  spectral  peak  at  the  first 
discrete  frequency  in  the  calculated  spectrum  (bin  1) . 

The  sampled  complex  return  signal  after  removal  of 
the  cable  profile,  as  before,  is 
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where 


d  =  phase  difference  between  the  two  channels; 
ideally,  0-  =  90  degrees 
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In  addition,  if  there  is  an  amplitude  imbalance  between  channels 
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where 


a  =  magnitude  of  the  real  (in-phase)  part  of  the 
return  signal 


0  =  magnitude  of  the  imaginary  (quadrature)  part  of 
the  return  signal,  ideally,  a  =  0. 


The  transformation  to  frequency  space  yields: 
i  N-l 
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*V  [a  cos  ^  +  j8  cos(^S  +  *)] 
n=0  N  N 


.  .  2irn 
exp  (~D  -jj-  k) 


(5.9) 


Substituting  the  realtion  cos  <J>  «  ^ |^exp ( j <j> )  +  exp(-j<J>)J  yields 

,  N-l  N-l  ,  N-l  _  „ 

F  |x (n)  |  =  |  l  exP(-j  ^(k-1))  +  f  £  exp (- j^(k+l) ) 

1  ’ n=0  *  n=0  N  1  n=0  n 


N-l 


+  ^  exp(je)  l  exp (- j  (k-1) ) 
z  n*0  N 


N-l 


+  ^  exp(-je)  l  exp (- j  (k+1) ) 

z  n=0  N 
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%  6(k-l)  +  ^  6(k+l)  +  ^  exp(je)  o(k-l) 


+  ^  exp(-je)  6 (k+1) 


(5.10) 


Gathering  terms  gives; 


F  (x(n)}N  =  6(k-l)[|  +  #  exp(je)] 

1  J  n=0 


+  6 (k+1) ^  exp(-je)] 


j  <S(k-l)[a  -  3  sin0  +  jB  cos0] 


+  ^  6 (k+1) [a  +  B  sin0  +  jB  cos0]  (5.11) 


Let  y  =  ratio  of  magnitudes  of  the  positive  and  negative 
frequency  components 

|  a  -  6  sinO  +  jB  cos  0| 

| a  +  B  sine  +  j  cos  0 | 


or,  in  dB, 


Y 


10  log10 


+  B  -  2aB  sin  0 

__  j 

a  +  3  +  2a3  sxn  0 


(5.12) 


In  the  Queen's  guided  radar  system  the  measured 
errors  of  amplitude  (1  dB)  and  phase  (3.0  degrees)  give  a 
Y  of  -27  dB.  The  approximate  value  of  y  seen  in  the  signatures 
in  figures  5.3  through  5.11  is  -25  dB. 


5.6  Summary 

It  has  been  demonstrated  that  the  target  return  signal 
contains  information  about  the  kinematics  of  motion  of  the 
target.  This  information  could  be  useful  for  classifying 
targets.  It  is  contained  in  the  portions  of  the  return  signal 


which  are  doppler  shifted  due  to  the  instantaneous  velocity 
of  each  physical  section  of  the  target. 

The  appropriately  filtered  and  windowed  fourier 
transform  of  the  return  signal  is  referred  to  as  the  target 
signature.  In  all  signatures  a  large  peak  was  evident  which 
corresponded  to  the  average  velocity  of  the  target.  In 
addition ,  secondary  peaks  were  observed  due  to  motion  of 
the  legs  or  arms  of  a  human  target.  The  position  of  these 
side  lobes,  with  respect  to  the  main  peak  representing  average 
velocity,  reflects  characteristics  of  normal  human  locomotion. 
The  target  signature  of  a  dog  is  quite  different  from  that  of 
a  human  because  of  the  different  pattern  of  motion  of  the 
four  legs  of  the  dog. 

A  simple  target  such  as  an  automobile  or  the  wooden 
trolley  produced  a  signature  which  had  no  noticeable  features 
other  than  the  principal  peak. 

Some  inaccuracies  in  the  receiver  were  observed  to 
introduce  spectral  leakage  to  the  image  frequency  giving  rise 
to  anomalous  peaks  in  the  target  signature.  However,  these 
inaccuracies  could  be  compensated  for  with  slightly  more 
sophisticated  signal  processing  techniques. 

Much  work  has  been  done  on  target  signature  analysis 
for  other  types  of  detection  systems.  The  signal  pro¬ 
cessing  techniques  that  have  been  developed  could  well 
be  applied  to  the  guided  radar  case.  However,  the  signal 
processing  would  have  to  be  extended  to  deal  with  signifi¬ 
cant  differences  between  the  field  configuration  in  normal 
radar  and  that  in  guided  radar. 

The  return  signal  and  therefore  the  target  signature 
is  a  result  of  the  interaction  of  the  complex  target  with 
the  electromagnetic  field  produced  by  the  radar  system.  In 
a  normal  radar  this  field  is  spatially  uniform.  In  guided 
radar  the  target  is  interacting  with  a  distinctly  non- 
uniform  field  distribution.  The  contribution  of  each 
section  of  the  target  to  the  total  signature  is  highly  de¬ 
pendent  upon  its  exact  position  with  respect  to  the  sensors. 
More  experimental  and  theoretical  analysis  would  be  required 
to  determine  the  appropriate  techniques  to  deal  with  the 
added  complications  of  the  non-uniform  field.  One  approach 
would  be  to  simulate  the  target  return  signal  and  the  target 
signature  based  upon  known  field  configuration  data  and 
the  kinematics  of  motion  of  the  target.  Initial  experiments 
have  been  performed  using  this  approach  and  are  discussed 
in  the  next  section. 


6  Simulation  of  Target  Signatures 


6.1  Introduction 

An  attempt  was  made  to  reproduce  the  form  of  the 
target  signatures  obtained  experimentally  for  human  targets 
by  modelling  the  human  as  a  collection  of  targets  all  moving 
with  the  same  average  velocity  but  each  having  a  different 
cross  section  and  instantaneous  velocity.  The  model  was 
intended  to  demonstrate  the  characteristics  of  human  loco¬ 
motion  which  produce  the  distinctive  features  in  the  signature 
described  in  the  previous  chapter. 


6.2  Return  Signal,  single  component  target 


Initially,  the  target  was  assumed  to  be  a  single 
point  object  moving  uniformly  along  the  cables  in  the  x 
direction.  The  orientation  of  the  spatial  axes  with  respect 
to  the  leaky  cable  sensor  is  shown  in  figure  6.1  and  is: 


x  axis:  parallel  to  the  cables 

y  axis:  perpendicular  to  the  cables 

z  axis:  perpendicular  to  the  plane  of  the  cables 


The  point  target  is  assumed  to  have  a  reflectivity  o  and  to 
be  at  position  (xQ,y0,zo)  at  time  t  =  0.  At  any  other  time 
t,  the  target  position  is  (x,y,z)  and  its  instantaneous 
velocity  is  v(t).  If  vx,  Vy  and  vz  are  the  x,  y,  and  z 
components  of  v,  respectively,  the  equations  of  motion  of 
the  point  object  can  be  written: 


x(t) 


t 


v  dt 
x 


y  (t) 


v  dt 

y 


z  (t) 


v  dt 
z 


(6.1) 

(6.2) 


(6.3) 


Both  the  I  and  Q  components  of  the  target  vector  are 

2tt f  x 

a  sinusoidal  function  of  - 

v 

coax 

where:  f  =  carrier  frequency 

Vcoax  =  velocity  of  propagation  of  the  carrier  in  the 
coaxial  cable. 
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In  addition, each  component  will  be  amplitude  modulated  by 
a  function  of  position  which  represents  the  varying  sensi¬ 
tivity  of  the  sensor  in  the  three  spatial  directions.  For 
convenience  this  function  will  be  divided  into  two  compo¬ 
nents:  one,  S(x)  representing  the  sensitivity  variation 
along  the  length  of  the  cables,  and  the  other,  S(y,z), 
representing  the  sensitivity  in  the  other  two  dimensions. 

From  field  site  experiments  it  is  known  that  S(x) 
is,  in  general,  a  slowly  varying  function  of  x.  The  magni¬ 
tude  of  S(x)  depends  upon  the  position  of  the  cables  with 
respect  to  the  ground  as  does  the  spatial  frequency.  The 
variation  of  S(x)  would  not  be  expected  to  greatly  influence 
the  target  signature  since  it  would  transform  to  a  low 
value  in  frequency  space.  S(x)  will  be  assumed  to  be  a 
constant. 

The  sensitivity  in  the  plane  perpendicular  to  the 
length  of  the  cables  is  distinctly  non-uniform.  The  res¬ 
ponse  is  determined  by  factors  such  as  the  distance  between 
the  cables,  the  type  of  cable  and  its  position  relative  to 
the  ground,  and  the  frequency  of  operation  of  the  system. 

The  sensitivity  contours  in  the  y-z  plane  have  been  measured 
[  11  ]  and  have  also  been  roughly  estimated  from  the  theoretical 
power  density  distribution  for  a  2  wire  d.c.  transmission 
line.  Initially,  the  sensitivity  in  the  transverse  plane  will 
be  assumed  to  be  constant.  Thus  the  model  will  not  include 
effects  due  to  the  lifting  of  the  feet  during  each  step  or 
due  to  deviation  of  the  target  in  the  y  direction  from  the 
midline  between  the  cables. 

The  complete  analog  baseband  signal  produced  by  a 
single  point  object  is: 

2irf 

I(t)  =  a  S(x)S(y,z)  sin(- - —  x)  (6.4) 

coax 


2TTf 

Q  ( t)  =  o  S  (x)  S  (y ,  z)  cos  (— - —  x)  (6.5) 

coax 


Multiple  Component  Targets 


Consider  a  target  consisting  of  n  distinct  com¬ 
ponents  having  target  cross-sections  {o-^,  o 2,  ...»  on} 
respectively.  In  addition  let  the  instantaneous  velocity 
of  these  components  at  time  t  be  (v^,  V2,  .  .  .  ,  vn)  res¬ 
pectively.  Since  all  n  components  are  part  of  one  composite 
body,  they  all  have  the  same  average  velocity.  Each  com¬ 
ponent  will  be  considered  to  be  a  point  object.  At  time 
t  =  0  let  the  position  of  the  nth  component  be  (xno,  yno,  zn„) . 
Equations  6.4  and  6.5  can  now  be  used  to  model  the  contribution 
of  each  component  of  the  target  to  the  total  return  signal. 
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I 

I 


I(t)  =  a,  S  (x, )  Sty^)  sin  - -  x 

A  x  coax 

2irf 

+  On  S  (x«)  S  (y_  ,  z5)  sm 

^  z  z  z  vcoax 


+  °n  S(V  Slyn'zn> 


sin 


2irf 
_ c_ 

v 

coax 


2*f 

Q(t)  =  o,  S (x, )  S(y1,z1)  cos  - - — 

1  A  A  coax 


2irf 

*  °2  S(W  ^ 


+  on  S(xj)  S(yn,*n> 


2lTf 


cos 


coax 


Hence 


n  2rri 

Kt)  -  I  o.  S(x  )  S(y.,z.)  sin  — 

i=l  X  CO 


n  2irj 

Q(t)  =  I  o.  S(x.)  S{yi,zi)  cos  — 
i-1  co 


where 


z . 
1 


(6.12) 


dt 


Equations  6.8  through  6.12  form  the  complete  simu¬ 
lation  of  an  ideal  return  signal  from  a  target  with  known 
motion-characteristics.  In  a  practical  system  this  basic 
return  signal,  along  with  extraneous  noise  caused  by  terrain 
and  weather  conditions,  and  system  noise,  will  determine  the 
actual  return  signal. 

Once  the  baseband  return  signal  is  generated  on  the 
basis  of  the  above  model  and  the  kinematics  of  motion  of 
the  target,  it  is  processed  in  a  fashion  identical  to  that 
used  to  obtain  a  target  signature  from  an  experimentally 
obtained  return  signal  (see  figure  5.1). 


6.4  Results 

6.4.1  Simple  Target 

Figure  6.2(a)  shows  the  signature  for  a  single 
component  simulated  target  moving  at  a  uniform  velocity 
of  0.7  m/sec.  This  would  represent  the  signature  of  a 
human  target  when  the  trunk,  moving  at  the  average  velocity, 
is  the  major  body  part  causing  the  target  vector.  It  would 
correspond  to  the  situation  when  the  leaky  cables  are  raised 
above  the  ground  so  that  the  motion  of  the  legs  would  have 
little  or  no  influence  on  the  return  signal  and  when  the 
arms  are  held  stationary  with  respect  to  the  body.  For 
comparison  figure  6.2(b)  shows  the  signature  for  a  human 
target  when  the  cables  were  at  a  height  of  2  m,  well  away 
from  the  leg  motion.  In  both  cases  there  is  a  large  peak  in 
the  signature  corresponding  to  the  average  target  velocity  but 
little  other  distinct  structure. 


6.4.2  Three  Component  Target 

To  include  the  effects  of  the  motion  of  the  legs  in 
the  simulated  target  signature  a  three  component  target  was 
used.  One  component  of  the  target  moved  at  a  uniform  velocity 
and  represented  the  trunk  and  upper  limbs  of  the  human.  Each 
of  the  other  two  components  represented  the  lumped  effect  of 
one  leg  and  foot.  Naturally,  they  moved  at  the  same  average 
velocity  as  the  first  component;  however,  they  had  an 
instantaneous  velocity  which  varied  with  time  in  the  same 
manner  as  the  instantaneous  velocity  of  the  human  heel 
during  normal  walking.  The  legs  were  assumed  to  be  180 
degrees  out  of  phase  with  each  other ^ 

The  velocity  time  curves  used  in  the  model  to  represent 
the  leg  motion  are  shown  in  figure  6.3.  They  were  obtained  in 
a  study  by  Winter  [  12 ]  of  the  kinematics  of  normal  human  loco¬ 
motion.  The  curves,  in  fact,  are  for  the  motion  of  the  heel  of 
the  foot. 
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Figure  6.2:  Target  Signatures  -  Human  Target 
Cables  in  Air 


A  signature  for  a  walking  human  target,  modelled  as 
three  lumped  components,  is  shown  in  figure  6.4(a).  Figure 
6.4(b)  is  a  corresponding  signature  experimentally  obtained 
with  the  cables  on  the  ground.  A  comparison  of  the  two 
signatures  reveals  that  the  simulated  target  signature  shows 
some  features  characteristic  of  the  experimentally  obtained 
signature.  In  particular,  the  distance  in  frequency  space 
between  adjacent  peaks  is  approximately  three  times  the 
frequency  resulting  from  the  average  velocity  in  both  sig¬ 
natures  . 

To  produce  the  simulated  signature  the  reflection 
coefficient  (or  cross-section  6)  of  the  legs  was  arbitrarily 
chosen  to  be  1/10  that  of  the  body  and  no  attempt  was  made 
to  match  the  peak  heights  in  the  simulated  signature  to  those 
in  the  experimentally  obtained  one.  Also  the  system  inaccur¬ 
acies  discussed  in  section  5.5  which  cause  extra  peaks  to  occur 
in  the  experimental  signature  were  not  included  in  the  simul¬ 
ation. 


6.5  Conclusions 

A  very  simplified  model  of  the  human  body  as  a  target 
gives  signatures  which  correspond  in  gross  detail  to  those 
obtained  experimentally.  The  use  of  the  motion  characteristics 
of  the  feet  as  representative  of  those  of  the  legs  may  explain 
some  of  the  difference  between  the  synthesized  and  experimentally 
obtained  signatures.  A  more  exact  model  would  incorporate  target 
components  (with  appropriate  motion  characteristics)  corresponding 
to  the  different  sections  of  the  leg.  Details  of  the  field  dist¬ 
ribution  (S(y,z))  and  target  component  motion  in  the  y  and  z 
directions  should  also  be  included  in  a  complete  simulation  of 
a  target  signature. 

These  initial  results  indicate  that  it  should  be 
possible,  in  principle,  to  create  a  library  of  target  signatures 
on  the  basis  of  which  general  target  classification  could  be 
performed.  A  knowledge  of  the  dynamics  of  target  motion  could 
be  used  to  predict  distinctive  features  of  expected  signatures. 

The  experiments  in  which  the  signature  for  human  targets  was 
measured,  reported  in  Section  5,  were  conducted  in  a  variety 
of  different  weather  conditions  including  both  clear  dry 
weather  and  heavy  rain.  Even  though  the  system  response  changes 
under  different  environmental  conditions,  the  changes  are  such 
that  they  do  not  greatly  affect  major  features  of  the  signatures. 

Naturally,  humans  can  dramatically  change  the  dynamics 
of  their  locomotion  (walking  upright,  crawling  etc) .  Animals 
on  the  other  hand  are  not  likely  to  do  so.  Therefore,  the 
recommended  approach  to  reducing  false  alarms  caused  by  animals 
would  be  to  identify  those  characteristics  of  motion  of  animals 
which  are  distinct  and  relatively  constant  and  which  result  in 
clear  patterns  in  the  target  signatures.  A  signature  simu¬ 
lation  technique  could  be  useful  in  identifying  these  distinc¬ 
tive  signature  patterns. 
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Signature  simulation  could  also  be  a  useful  method 
for  assessing  the  effective  cross-section  of  various  parts 
of  the  human  body.  This  would  require  a  knowledge  of  the 
field  distribution  6 (x)  and  6(y,z). 

Although  comments  in  this  report  have  been  mainly 
limited  to  longitudinal  motion  of  the  target  distinct  target 
signatures  should  also  be  produced  by  transverse  crossings 
of  the  sensors.  The  characteristics  of  these  signatures 
would  result  from  the  interaction  of  the  target  with  the 
changing  field  strength  in  the  y  direction  rather  than  from 
an  actual  doppler  frequency  shift  in  the  return  signal. 


7  Single  Cable  System 


7 . 1  Introduction 

In  all  the  work  discussed  so  far  in  this  report  two 
leaky  cables  have  been  used  as  the  sensor  in  a  guided  radar 
system  and  targets  have  been  detected,  located  or  identified 
on  the  basis  of  the  interaction  between  the  target  and  the 
electromagnetic  field  between  the  cables.  This  configuration 
is  well  suited  to  perimeter  surveillance  since  the  detection 
zone  is  confined  to  the  area  around  the  cables.  However,  in 
principle,  it  should  be  possible  to  perform  the  same  functions 
with  a  centrally  located  antenna  and  a  single  surrounding  leaky 
cable.  In  effect,  a  distributed  antenna  is  replaced  with  a 
point  source  radiator.  The  wave  propagating  from  the  central 
antenna  couples  into  the  leaky  cable.  A  target  in  the 
vicinity  of  the  cable  or  between  the  cable  and  the  antenna 
influences  the  amount  of  coupling  and  therefore  changes  the 
signal  monitored  at  one  end  of  the  cable.  Through  this  change 
the  target  can  be  detected. 

In  a  two  cable  system  the  location  of  the  target  can 
be  determined  by  a  technique  similar  to  time  delay  reflecto- 
metry.  R.F.  energy  is  transmitted  along  the  sensor  in 
pulses.  The  target  location  is  determined  by  the  time  delay 
between  the  pulse  transmission  and  the  appearance  of  the  target 
vector.  In  the  single  cable  short  perimeter  system  a  similar 
technique  could  be  used  with  suitable  modification  to  account 
for  the  fact  that  a  spherical  wavefront  is  being  transmitted 
from  a  central  antenna.  Transmission  of  a  simple  continuous 
wave  simplifies  the  processing  system  but  eliminates  the 
location  capability  for  a  perimeter  of  useful  length.  However, 
by  a  simple  modulation  technique  the  target  location  capability 
can  be  included  in  a  CW  system. 

The  purpose  of  the  experiments  reported  in  this  section 
was  to  investigate  the  performance  of  leaky  coaxial  cables  in  a 
single  cable  system.  Furthermore,  a  technique  was  investigated 
which  would  allow  the  location  of  the  target  to  be  determined 
even  when  operating  in  a  CW  mode. 

Section  7.2  gives  a  brief  analysis  of  the  operation 
of  the  system  and  contains  a  discussion  of  the  constraints  on 
the  proposed  location  technique.  The  hardware  and  software 
are  described  in  Section  7.3.  The  experiments  performed  and 
the  results  are  presented  in  Section  7.4.  Section  7.5 
presents  a  discussion  of  the  effects  of  errors  in  the  analog 
hardware  on  the  accuracy  of  the  system  and  Section  7.6  is  a 
summary . 


7.2  System  Analysis 

7.2.1  Basic  Operation 

The  system  to  be  analysed  consists  of  an  antenna  placed 
in  the  centre  of  a  perfect  circle  of  leaky  coaxial  cable.  The 
cable  has  a  receiver  connected  to  one  end  and  a  matched  load 
connected  to  the  other  end.  The  antenna  is  connected  to  a 
transmitter  which  is  producing  a  signal  at  a  frequency  fQ. 

The  setup  is  shown  in  figure  7.1. 
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Figure  7.1:  System  Configuration 


The  electromagnetic  field  produced  by  the  antenna  is 
assumed  to  be  symmetric  with  respect  to  azimuthal  angle. 

Since  the  distance  from  the  antenna  to  the  cable  is  the 
same  for  all  points  on  the  cable,  the  phase  of  the  signal 
at  the  cable  is  also  the  same  for  all  points  on  the  cable. 

A  portion  of  the  electromagnetic  field  at  the  cable  is  con¬ 
verted  into  a  coaxial  mode  and  propagates  down  the  leaky 
cable  to  the  receiver  or  to  the  load. 

If  a  target  enters  the  detection  zone  around  the 
cable  the  field  will  change  as  will  the  amount  of  energy 
coupled  into  the  cable.  The  change  in  voltage  at  the  receiver 
can  be  shown  to  be:  (see  Appendix  B) 


V  =  AM  sin(2irf  t  +  0  +  $  )  7.1 

T  x  o  x 

where 

f  is  the  carrier  frequency 

6  is  a  phase  delay  due  to  the  propagation  time  from 
the  antenna  to  the  cable.  It  is  constant  for  all 
points  on  the  cable 

@  is  a  phase  delay  due  to  the  propagation  time  along 
x  the  coaxial  cable  from  the  targets'  location  to 
the  receiver 

AM  is  the  change  in  the  signal  coupled  into  the  cable 
at  the  location  of  the  target 


The  original  transmitted  signal  was  assumed  to  be  of  the  form 
sin(2TrfQt)  . 

The  presence  of  a  target  is  indicated  by  AMX  ^  0 

The  phase  term  3X  is  directly  proportional  to  the  distance 
along  the  cable  from  the  target  to  the  receiver 


6 


x 


27TXf0 


V 

coax 


7.2 


where  x  is  the  distance  along  the  cable  from  the  target  to 
the  receiver 

vcoax  is  the  velocitY  propagation  in  the  cable 


The  length  of  the  cable  over  which  unambiguous  location 
is  possible  is  restricted.  This  restriction  is  discussed  in  the 
next  section. 
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7.2.2  Cable  Length  Restriction 

The  form  of  equations  7.1  and  7.2  places  a  severe 
restriction  on  allowable  perimeter  length  for  unambiguous 
target  location.  The  phase  term  @x  must  be  less  than  tt 
radians  since 

AM^  can  be  positive  or  negative  and  >s 


sin(2irf  t  +  0  +  6)  =  -  sin(27rf  t  +  0  +  8  +  n) 

ox  ox 


Thus 


v  *  coax 
x  -  ~2 1 

o 


X  =  wavelength  of  the  electromagnetic  field  in  the 
coaxial  cable 

For  example,  for  a  perimeter  of  125  m  and  vcoax  =  0.8c 

f  must  be  less  than  960  kHz.  Standard  radiax  leaky  cables 
are  very  inefficient  at  this  low  frequency. 

7.2.3  Modulation  Technique 

A  more  suitable  operating  frequency  for  Radiax  is  in 
the  VHF  range  (50  MHz  to  150  MHz) .  If  a  carrier  frequency 
in  this  range  is  chosen, the  ability  to  locate  targets  can  be 
maintained  by  modulating  the  carrier  with  a  frequency  chosen 
on  the  basis  of  the  required  perimeter  length. 

Double  sideband-suppressed  carrier  modulation  was 
chosen.  The  transmitted  waveform  is  now 


V  =  sin(2irf  t)sin(2irf  t) 
o  m 


where 


fQ  is  the  carrier  frequency 
fm  is  the  modulation  frequency 


It  can  be  shown  by  a  simple  extension  of  the  analysis  in 
Appendix  B  that  the  target  vector  now  has  the  form 

VT  =  AMxsin<2,lfC)t  +  0  +  8x)sin(27rfmt  +  <P  +  Yx>  ’ 

The  phase  terms  6  and  $  are  constants  determined  by  the  geo¬ 
metry  of  the  system  and  do  not  influence  the  analysis.  They 
will  be  ignored  from  this  point  on.  Synchronous  detection 
with  appropriate  filtering  at  both  the  carrier  frequency  and 
the  modulation  frequency  yields  four  signals  containing  the 
desired  target  information 
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AMsin8xsinvx 
AMsingxcosYx 
AMcosgxsinyx 
AMcosg  cosy 

X  X 


Target  vector  magnitude  is  given  by 


mt  -  /  AI J  +  AQ-  +  AI-  +  AQ2 


Target  vector  phase  can  be  calculated  from 


Y  =  tan 
x 


-1  AQ1 


sinB  ^  0 
x 


cos&  t  0 
x 


The  carrier  frequency  can  be  chosen  high  enough  for 
efficient  coupling  into  the  coaxial  cable  and  the  modulation 
frequency  low  enough  to  provide  unambiguous  target  location 
over  the  desired  perimeter. 


7.3  System  Description 

7.3.1  Hardware 

A  block  diagram  of  the  hardware  used  in  the  single 
cable  system  is  shown  in  figure  7.2.  The  transmitter  con¬ 
sisted  of  a  60  MHz  carrier  oscillator  amplitude  modulated 
by  a  700  kHz  modulation  oscillator,  and  a  five  watt  power 
amplifier . 

The  receiver  consisted  of  two  stages  of  synchronous 
detection  each  followed  by  the  appropriate  filtering.  The 
first  stage  detector  synchronously  demodulated  the  60  MHz 
carrier  frequency.  A  hybrid  was  used  to  provide  the  90° 
phase  shift  of  the  reference  carrier  signal  for  the  quadra¬ 
ture  mixer.  The  resulting  in  phase  (I)  and  quadrature  (Q) 
components  from  the  first  detector  were  filtered  through  a 
pair  of  9  pole  low-pass  filters  with  a  cut-off  frequency  of 
700  kHz. 

The  I  and  Q  components  from  the  first  detector  were 
then  separately  processed  by  a  pair  of  synchronous  detectors 
using  the  700  kHz  modulation  oscillator  as  the  reference 
signal.  In  this  case,  specially  cut  delay  lines  were  used 


Figure 


to  shift  the  reference  oscillator  by  90° .  The  four  resulting 
signals  were  filtered  by  low  pass  filters  having  a  cut-off 
frequency  of  30  Hz  and  then  sequentially  sampled  and  con¬ 
verted  to  a  digital  representation.  All  four  signals  were 
digitized  in  16  ms  giving  an  effective  profile  sampling  rate 
of  60  Hz. 

7.3.2  Software 

A  flow  diagram  for  the  processing  scheme  is  shown  in 
figure  7.3.  Each  time  sample  of  the  return  signal  consisted 
of  four  numbers  -  digitized  versions  of  1^,  Q^,  I2  and  Q2» 

The  initial  value  of  each  of  these  signals  was  used  as  the 
reference  profile  -  the  system  response  with  no  target  present 
within  the  detection  zone  of  the  system. 

The  reference  profile  was  subtracted  from  each  sub¬ 
sequent  sample  of  the  profile  to  form  a  target  response  vector 
(AI^,  AQ^ ,  Al2'  ’  T*ie  ma9nitu<*e  the  target  response 

vector  was  then  calculated  M  =  /ai^  +  AQ^  +  AI^  +  AQ^  . 

If  this  magnitude  was  below  a  threshold  level,  chosen  on  the 
basis  of  system  noise  levels,  it  was  assumed  that  no  target 
was  within  the  detection  zone  and  the  next  profile  sample 
was  processed.  If  the  target  vector  magnitude  was  greater 
than  the  threshold  a  comparison  was  made  between  and  M2 

where  =  / AI ^  +  AI2  M2  =  /AI2  +  AQ2  .  The  location  of 

the  target  was  then  calculated  on  the  basis  of  a  target  phase 
-1  AQn 

Yx  =  tan  where  n  =  1,2  depending  upon  which  of  or 

n 

2  was  largest. 


7.4  Results 

7.4.1  Standard  Walk 

Most  of  the  measurements  on  the  system  were  made  by 
means  of  a  standard  walk  in  which  a  human  target  followed  a 
fixed  path  around  the  protected  perimeter.  This  path  is 
shown  in  figure  7.4. 

At  the  start  of  the  experiment  the  target  was  well  out 
side  the  detection  zone  assuring  a  reliable  recording  of  the 
reference  profile.  The  target  crossed  the  cable  at  one  end 
and  walked  at  a  constant  velocity  approximately  35  cm  from 
the  cable  on  the  inside  of  the  loop. 


7.4.2  Experiments 

A  variety  of  experiments  were  performed  with  different 
system  configurations.  These  are  summarized  in  Table  7.1 
Each  experiment  was  performed  a  number  of  times  and  typical 
results  will  be  discussed  below. 


Figure  7.3:  Processing  Flow  Diagram 
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Standard  Walk  Path 
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Table  7 . 1 


SUMMARY  OF  SINGLE  CABLE  EXPERIMENTS 


Cable 

Type 

Transmitter 

Cable 

Position 

Special 

Coriditions 

RX4-3 

Antenna 

ground 

Cable 

ground 

Antenna 

1  metre 

Antenna 

ground 

large  stationary 
target 

Antenna 

ground 

radial  walk 

RX4-1 

Antenna 

ground 

Cable 

ground 
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7.4.3  Basic  System  Performance 

In  figure  7.5  are  plotted  the  results  of  a  typical 
standard  walk  experiment.  An  RX4-3  cable  was  used  on  the 
ground.  The  lower  curve  is  a  plot  of  the  target  vector  mag¬ 
nitude  Aw  versus  time.  Vne  upper  curve  is  a  plot  of  the 
calculated  target  positon  versus  time.  The  straight  line 
is  the  best  fit  straight  line  to  the  data  using  a  minimum 
root-mean-square  error  criterion.  Marked  on  the  graph  is 
the  peak  value  of  the  target  response.  The  results  plotted 
in  figure  7.6  are  for  a  similar  experiment  with  the  exception 
that  an  RX4-1  cable  was  used. 

It  is  clear  that  the  scheme  for  target  location  works 
well.  The  ripple  on  the  target  position  plots  results  from 
differences  in  the  gain  of  the  amplifiers  in  the  four  channels 
and  deviation  from  perfect  quadrature  in  the  synchronous 
detectors.  This  will  be  discussed  further  in  section  7.5. 

The  target  vector  magnitude  plots  are  similar  in  detail  with 
the  exception  that  the  target  response  for  RX4-1  is  15  dB  lower 
than  that  for  RX4-3. 

7.4.4  Performance  with  Cable  as  Transmitter 

The  reciprocity  theorem  implies  that  it  should  be 
possible  to  reverse  the  positions  of  the  cable  and  the  central 
antenna.  Figure  7.7  is  a  plot  of  the  results  from  a  standard 
walk  with  the  transmitter  connected  to  the  cable  and  the 
receiver  connected  to  the  central  antenna.  The  performance 
is  essentially  the  same  as  when  the  central  antenna  was  used 
as  the  transmitter. 


7.4.5  Elevated  Cables 

For  a  number  of  walks  the  cable  was  suspended  by 
plastic  clips  on  wooden  stakes  approximately  1  metre  above 
the  ground.  Figure  7.8  shows  the  results  of  one  such  experi¬ 
ment.  The  most  obvious  change  is  in  the  target  magnitude 
plot.  Large,  regularly  spaced  nulls  appear  where  the  system 
sensitivity  becomes  very  low.  At  these  nulls  the  target 
location  response  also  becomes  inaccurate. 

The  sensitivity  nulls  most  likely  occur  because  of 
interference  between  several  modes  of  propagation  along 
the  cable.  One  mode  would  be  a  coaxial  mode  within  the 
cable  and  a  second  mode  would  propagate  between  the  outer 
conductor  of  the  cable  and  the  ground  plane  produced  by  the 
earth.  These  modes,  propagating  at  different  velocities, 
interfere  to  produce  the  senstivity  nulls.  Similar  periodic 
sensitivity  nulls  have  been  observed  in  a  normal  two  cable 
guided  radar  system  when  the  cables  were  suspended  in  the 
air  r 13,14]  and  a  theoretical,  coupled-mode  model  has  been 
used  to  determine  the  parameters  which  influence  the  spacing 
and  depth  of  the  nulls  in  a  two  cable  system  [15]. 

The  decrease  in  location  accuracy  at  the  nulls  is  a 
result  of  the  decreased  signal- to-noise  ratio  at  these  points. 
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Target  Response  (volts)  Target  Position  (m) 


Target  Response  (volts)  Target  position  (m) 


7.4.6  Performance  with  Large  Stationary  Obstacles 

A  series  of  experiments  was  performed  to  determine 
the  influence  of  a  large,  stationary,  metallic  object  (an 
automobile)  on  the  performance  of  the  system.  The  car  was 
parked  inside  the  loop  of  cable  with  the  side  of  the  car 
approximately  1-1/2  feet  from  the  cable.  The  target’s  (a 
human)  path  was  the  same  as  in  the  standard  walks  with  the 
target  passing  between  the  car  and  the  cable. 

A  typical  result  is  shown  in  figure  7.9.  The  overall 
system  sensitivity  was  about  the  same  as  it  would  be  with  no 
obstacle  present  except  when  the  target  was  "shadowed"  by  the 
car.  At  that  point  the  sensitivity  dropped  by  about  3  dB.  The 
target  location  accuracy  was  not  affected. 


7.4.7  Radial  Walks 

The  last  experimental  series  involved  a  radial  walk. 

The  target  started  approximately  10  m  outside  the  loop  and 
walked  directly  towards  the  antenna.  The  target  remained 
at  the  antenna  for  fifteen  seconds,  turned,  and  followed 
the  same  radial  path  back  to  the  starting  point.  Figure  7.10 
is  a  plot  of  typical  data  collected  during  such  a  radial  walk. 

The  target  is  clearly  visible  in  both  the  target  res¬ 
ponse  plot  and  the  target  location  plot.  As  the  target 
approached  the  antenna  the  sensitivity  dropped  below  the 
threshold  for  detection  chosen  for  this  experiment  (-8  dB) . 

When  close  to  the  antenna  the  target  again  was  detected 
although  it  was  not  accurately  located.  A  very  similar 
result  was  produced  when  the  cable  was  connected  to  the 
transmitter . 

7.5  System  Simulation 

7.5.1  Introduction 

On  the  basis  of  the  analysis  of  the  system  (section  7.3 
and  Appendix  B)  it  is  possible  to  simulate  the  return  signal 
produced  at  the  receiver  when  a  target  is  within  the  detection 
zone  of  the  system.  Furthermore,  the  processing  that  occurs 
in  the  analogue  portion  of  the  receiver  (see  figure  7.3)  can 
also  be  simulated.  If  the  resulting  digital  representation 
of  the  signal  is  then  processed  by  the  normal  target  detection 
and  location  algorithm,  it  should  be  possible  to  determine  the 
effects  of  errors  in  the  analogue  processing  on  the  accuracy 
of  target  detection  and  location.  This  exercise  has  been 
carried  out. 

Under  ideal  conditions,  with  matched  gains  in  all 
four  analogue  channels  and  exact  90°  phase  shifts  in  the 
synchronous  detectors  the  algorithm  output,  with  a  simulated 
input  signal  is  shown  in  figure  7.11.  The  only  differences 
between  these  plots  and  the  ideal  straight  line  plots  results 
from  quantization  error  in  the  simulation  and  algorithm. 


7.5.2  Quadrature  Error 

A  phase  shift  of  other  than  90°  in  the  carrier 
detector  will  cause  a  change  in  two  of  the  four  signals 
produced  by  the  analogue  system.  Equations  7.5  can  be 
rewritten 


AI,  =  AMsin(6  +a)siny 
1  x  x 

AQ^  =  AMsin<ex+oi)  cosyx 

7.9 

AI~  =  AMcosg  siny 
z  x  x 

AQ_  =  AMcosg  cosy 

4  XX 


where  a  is  the  error  in  the  quadrature  oscillator. 

Since  target  location  is  calculated  either  from  AI./AQ^  or 
AI2/AQ2  the  location  plot  should  show  no  dependence  on1 
carrier  detector  quadrature  error.  The  error  will,  however, 
affect  the  target  magnitude  response.  Figure  7.12  illustrates 
these  points.  For  this  simulation  a  =  10°. 

An  error  in  the  baseband  detector  would  cause  equations 


7.5  to  have 

the  form 

AI1 

=  AMsin0x 

sin ( yx+6 ) 

AQi 

=  AMsin0x 

cosyx 

AI2 

=  AMcos0x 

sin(yx+6) 

AQ2 

=  AMcos  3 

X 

cosy 

V 

x 


where  6  is  the  error  in  the  baseband  detector. 

An  error  term  will  occur  in  both  the  target  magnitude  res¬ 
ponse  curve  and  the  target  position  curve.  The  form  of  the 
•error  is  illustrated  in  figure  7.13.  An  error  of  10°  was 
included  in  the  baseband  quadrature  oscillator. 


7.5.3  Amplifier  Gain  Imbalance 

Differences  in  the  gain  of  the  four  system  channels 
will  affect  both  the  target  magnitude  response  and  the  target 
location  accuracy.  The  precise  effect  of  an  error  will  depend 
upon  which  channels  are  miscalibrated . 
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then  the  magnitude  response  curve  will  be  affected  but  the 
target  location  plot  will  be  correct  since  target  location 
accuracy  depends  upon  AI^/AQ^. 

If  G  g  for  either  n  =  1  or  2 

n  yn 

then  both  the  magnitude  resoonse  and  the  location  accuracy 
will  be  influenced.  The  choice  between  calculating  target 
position  on  the  basis  of  AQ^/AIi  and  AQ2/AI2  depends  upon 
the  ratio  of  and  M?  mentioned  in  section  7.3.2.  As  the 
target  proceeds  around  the  cable  loop  the  algorithm  will 
alternate  between  using  the  first  pair  of  signals  and  the 
second  pair  of  signals  for  position  calculation.  In  a 
standard  walk,  therefore,  the  target  position  will  appear 
to  jump  back  and  forth  about  an  average  position  as  the 
algorithm  switches  from  using  the  first  pair  of  signals 
to  using  the  second  pair. 

Figure  7.14  shows  the  output  of  the  algorithm  when 
an  error  of  2  dB  has  been  introduced  into  the  gain  of  the 
AI^  channel. 

The  errors  introduced  into  the  simulated  analogue 
system  in  the  examples  shown  in  this  section  have  been 
severe  compared  to  what  can  be  achieved  in  practice.  From 
the  results  it  would  appear  that  neither  the  detection 
accuracy  nor  the  location  accuracy  of  the  system  is  parti¬ 
cularly  sensitive  to  precise  adjustment  of  system  parameters. 


7.6  Summary 

It  is  clear,  from  the  experiments  reported  in  this 
section,  that  it  is  possible  to  replace  one  of  the  leaky 
coaxial  cables  in  a  standard  guided  radar  system  with  a 
central  antenna  and  still  provide  intrusion  detection. 

When  the  cable  is  suspended  in  the  air  effects  are  seen 
which  are  similar  to  those  in  a  two  cable  system  due  to 
interference  between  several  modes  of  propagation  of  the 
energy  along  the  cable.  It  should  be  possible  to  minimize 
these  effects  by  the  same  techniques  as  have  been  suggested 
for  a  two  cable  system. 
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No  experiments  were  conducted  with  the  leaky  cable 
buried  in  the  ground.  However,  based  upon  experience  gained 
with  a  two  cable  system  both  with  the  cables  buried  and  on 
the  surface,  no  significant  change  in  the  performance 
would  be  expected  if  a  buried  cable  were  used. 

The  system  can  operate  with  either  the  cable  or  the 
antenna  as  the  transmitting  element.  The  target  detection 
performance  is  similar  in  both  cases.  However,  when  the 
cable  is  used  as  the  transmitting  element,  the  field 
strength  in  the  area  surrounding  the  protected  perimeter 
would  be  lower  since  a  large  portion  of  the  transmitter 
power  would  be  absorbed  by  the  matched  load.  This  would 
obviously  be  the  best  configuration  in  situations  where 
several  single  cable  detection  systems  must  operate  in 
close  proximity  to  each  other  or  where  radiation  from  the 
intrusion  detection  system  would  interfere  with  other 
systems . 

The  feasibility  of  the  proposed  modulation  scheme 
for  allowing  the  measurement  of  location  of  an  intrusion 
was  successfully  demonstrated.  The  maximum  perimeter  used 
was  125  metres,  not  because  of  any  inherent  restriction, 
but  simply  because  of  the  available  cable  lengths.  Peri¬ 
meters  of  any  circumference  could  be  protected  as  long  as 
the  modulation  wavelength  was  at  least  twice  the  circum¬ 
ference  and  the  signal  coupled  into  the  cable  was  large 
enough  to  be  detected.  The  accuracy  of  location  is  pro¬ 
portional  to  the  accuracy  of  measurement  of  a  phase  angle 
and  could  be  stated  as  being  the  appropriate  constant 
fraction  of  a  wavelength  of  the  modulation  signal.  A  higher 
modulation  frequency  therefore  gives  better  location  accuracy. 
The  best  choice  of  modulation  frequency  is  the  frequency 
having  a  wavelength  which  is  twice  the  length  of  the  protected 
perimeter.  Since  the  location  accuracy  depends  upon  modu¬ 
lation  wavelength,  as  long  as  a  good  signal  to  noise  ratio 
were  maintained  no  change  in  the  location  accuracy  would  be 
expected  if  the  cable  was  buried. 

Targets  could  be  detected  in  the  area  between  the 
cable  and  the  antenna.  Only  a  few  preliminary  experiments 
were  performed  to  investigate  the  target  detection  character¬ 
istics  within  the  loop.  It  appears  that  the  target  is  de¬ 
tected  within  an  annulus  about  the  cable  and  when  it  is  near 
the  central  antenna.  Target  location  is  accurate,  however, 
only  when  the  target  is  near  the  cable.  More  detailed 
experimentation  would  be  required  to  determine  whether  it 
would  be  possible  to  track  the  target  when  it  is  within  the 
protected  perimeter. 

The  simple  signal  processing  required  for  target 
detection  and  location  and  the  necessity  for  having  only  one 
cable  should  allow  the  development  of  a  practical  micro¬ 
processor  based  portable  short  perimeter  surveillance  system. 
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8 .  SUMMARY 


The  work  described  in  this  report  covers  a  number  of 
quite  different  aspects  of  the  performance  of  leaky  coaxial 
cables  in  a  variety  of  configurations.  Advances  have  been 
made  in  a  number  of  areas.  These  will  be  briefly  summarized 
be low . 


8.1  Environmental  Effects 

The  conditions  of  the  ground  around  the  leaky  coaxial 
cables  can  influence  their  performance  in  a  guided  radar 
system.  By  far  the  most  important  environmental  factor  is 
moisture,  either  in  the  form  of  rain  or  snow.  To  determine 
the  extent  of  the  influence,  long  term  monitoring  of  a  stable 
guided  radar  installation  must  be  undertaken.  The  results 
of  such  monitoring  have  been  reported  in  Section  3  of  this 
report.  The  winter  and  early  spring  seasons  were  chosen 
because  the  most  significant  changes  in  environmental  con¬ 
ditions  take  place  during  that  period. 

Sensitivity  changes  of  up  to  15  dB  were  noted  at  any 
point  on  the  sensor  although  the  average  range  of  sensitivity 
change  was  closer  to  10  dB.  The  presence  of  large  amounts 
of  moisture  decreased  the  sensitivity  while  small  amounts 
enhanced  sensitivity  by  1  or  2  dB.  The  change  in  sensitivity 
was  relatively  uniform  at  all  points  along  the  cable. 

Profile  varied  in  the  same  way  as  sensitivity  although 
to  a  greater  extent.  The  maximum  variation  for  any  cell  was 
22  dB  with  the  average  being  13  dB.  The  variations  of  pro¬ 
file  along  the  sensor  were  less  uniform  than  the  variations 
in  sensitivity. 

The  noise  in  the  system  was  relatively  constant 
throughout  the  experimental  period  and  was  generated,  in  the 
main,  in  the  receiver.  The  weather  conditions,  therefore, 
did  not  influence  the  system  noise  levels. 

Continued  long  term  monitoring  of  a  guided  radar 
installation  would  be  required  to  more  closely  delimit 
the  magnitude  of  the  environmental  effects  since  particular 
weather  conditions  occur  only  once  or  twice  during  a  season. 
This  monitoring  should  include  automated  measurement  of  soil 
temperature  and  relative  moisture  content.  Also  the  results 
presented  in  this  report  were  obtained  from  cables  laid  on 
the  surface  of  the  ground.  The  effects  of  weather  would  be 
expected  to  be  different  for  buried  cables.  For  example, 
a  greater  delay  should  occur  between  the  onset  of  precipi¬ 
tation  and  an  effect  being  noted  in  the  profile  or  sensiti¬ 
vity.  The  time  constants  of  the  effects  due  to  drying  of 
the  surface  soil  should  also  be  different.  A  facility,  such 
as  the  automated  data  logging  system  at  the  Queen's  Univer¬ 
sity  field  site  would  be  required  for  adequate  monitoring 
during  rapidly  changing  conditions. 

It  is  recommended  that  further  monitoring  of  an 
installed  leaky  coaxial  sensor  be  carried  out  to  achieve 
some  statistical  significance  in  the  results.  This  should 


be  done  for  a  variety  of  cable  installations  including  sur¬ 
face  mounted  and  buried  cables  and  perhaps  including  cables 
buried  in  soil,  gravel  and  sand. 

8.2  Velocity  Measurement 

The  results  presented  have  demonstrated  that  the 
velocity  of  a  target  moving  longitudinally  between  the 
sensors  can  be  measured,  based  upon  the  doppler  shift  in 
the  return  signal.  The  accuracy  of  the  measurement  depends 
upon  system  design.  The  velocities  of  two  targets  can  be 
measured  independently,  assuming  they  are  different  from 
each  other.  In  some  situations  different  signal  processing 
could  be  used  to  simplify  the  system.  For  example,  if  only 
one  target  is  present,  hard  limiting  and  zero-crossing 
counting  could  replace  the  fourier  transform  operation 
presently  used. 

8.3  Target  Classification 

The  experiments  reported  indicate  that  basic  target 
classification  can  be  performed  using  the  return  signal  in 
a  guided  radar  system.  The  classification  is  based  upon 
the  so-called  target  signature  -  the  log-magnitude  of  the 
fourier  transform  of  the  target  vector.  Only  rudimentary 
classification  has  been  demonstrated  and  a  large  amount  of 
work  in  this  area  remains  to  be  done.  The  results  of  much 
of  the  work  in  target  classification  in  normal  radar  sys¬ 
tems  can  be  applied  directly  to  the  leaky  cable  system 
although  there  are  peculiarities  of  the  guided  radar  system 
(non-uniform  field  distribution  around  the  sensor,  environ¬ 
mental  effects)  which  require  study  directly  within  the 
guided  radar  context  to  assess  their  influence  on  target 
signatures  and  target  classification. 

The  ability  to  classify  targets,  even  broadly,  has 
a  significant  effect  on  the  false  alarm  rate  of  a  deployed 
system.  It  would  also  allow  an  appropriate  response  to  a 
detected  intrusion.  For  these  reasons,  extension  of  the 
work  reported  here  is  highly  recommended. 

Only  classification  of  targets  moving  longitudinally 
between  the  sensors  was  attempted.  However,  because  of 
the  decidedly  non-uniform  field  distribution  in  the  trans¬ 
verse  plane  in  the  vicinity  of  the  cables  it  should  be 
possible  to  produce  a  "signature"  for  a  transverse  crossing 
and,  perhaps,  classify  on  the  basis  of  this  signature.  A 
study  of  this  possibility  is  also  recommended . 

8.4  Single  Cable  System 

The  performance  of  a  single  cable  system  with  a  centrally 
located  antenna  was  investigated.  Such  a  system  would  be 
satisfactory  for  protecting  relatively  short  perimeters. 

The  system  operated  in  a  CW  mode.  However,  a  modulation 
technique  was  included  which  allowed  the  position  of  the 
intrusion  to  be  measured. 


The  system  performed  well  and  its  feasibility  was  thus 
demonstrated.  The  area  used  for  the  experiments  was 
relatively  clear  of  stationary  obstacles  (except  when  an 
obstacle  was  intentionally  placed  near  the  cable) .  The 
influence  of  stationary  targets  should  be  studied  in  some 
detail . 

The  processing  involved  in  the  system  is  such  that 
it  is  possible  to  implement  it  on  a  microprocessor  based 
system  thus  resulting  in  a  compact  and  portable  signal 
processing  signal  processing  unit.  In  fact,  a  start  has 
been  made  in  this  direction  in  a  4th  year  project  performed 
last  year  at  Queen's  University.  (This  project  was  not 
funded  by  this  contact.)  An  algorithm  was  developed  for 
the  M6800  microprocessor  which  demonstrated  that  the 
processing  required  could  be  done  by  a  microprocessor.  An 
extension  of  this  work  should  be  contemplated. 
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APPENDIX  B 


Assume  that  the  coupling  from  the  external  field 
into  the  coaxial  cable  takes  place  at  only  three  discrete 
points  around  the  circumference  x,  y  and  z.  The  signal 
from  the  antenna  is  sin(2irf0t).  Three  voltage  components  will 
propagate  in  the  leaky  cable. 


Vx(t)  =  X  sin<2irfQt  +  6) 

Vy(t)  =  Y  sin ( 2nfQt  +  6)  B.l 

V  (t)  =  Z  sin(2nf  t  +  9) 
z  o 

where  f  is  the  carrier  frequency 

0  is  a  constant  caused  by  the  propagation  time  of 
the  signal  from  the  antenna  to  the  cable 

Xf Y , Z  are  "coupling  coefficients" 


The  coupling  coefficients  represent  the  proportion  of  the 
radiated  field  which  is  coupled  into  the  coaxial  cable  at  a 
particular  point.  In  an  ideal  system  with  uniform  sensitivity 
X,  Y  and  Z  would  be  equal. 

The  signal  at  the  receiver  is: 


V  (t)  =  X ' sin (2nf  t  +  0  +  0  ) 

K  OX 

+  y'sin(2trf  t  +  0  +  0,,) 
o  y 

+  Z'sin(2Tif  t  +  0  +  0  ) 
o  z 


B .  2 


where  X',  Y*  and  Z'  differ  from  X.  Y  and  Z  respectively,  due 
to  attenuation  in  the  cable 

0  ,  0  and  0  are  phase  terms  caused  by  the  propagation 
x  y  z 

delay  in  the  cable  from  the  particular  point  to  the 
receiver . 


Assume  that  a  target  enters  the  detection  zone  at  point  X. 
The  total  receiver  signal  will  now  be 


VRT 


=  (X 1 +AX* )  sin  ( 2irf  t  +  8  +  B  ) 

O  X 

+  Y  1  sin  (2trf  Qt  +  8  +  8y)  B.3 

+  Z 1  sin (2irf  t  +  8  +  B„) 
o  z 

The  change  in  received  signal  due  to  the  presence  of  the 
target  is 

VRT  '  VR  =  AX'  sin  (2TTfQt  +  8  +  8x>  B.4 

In  Section  7  VRT  -  VR  is  referred  to  as  the  target  vector  VT 
and  AX'  is  given  the  symbol  AM  . 
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MISSION 

of 

Rente  Air  Development  Center 

RAVC  plain  and  execute*  tie*  earch,  development,  tut  and 
selected  acquisition  programs  In  support  of  Command ,  Contact 
Communication*  and  Intelligence  (C*I)  activities.  Technical 
and  engineering  support  within  area*  of  technical  competence 
l*  provided  to  ESV  Program  Office*  (POs)  and  other  BSP 
elements.  The  principal  technical  mission  atiea*  axe 
communications,  electromagnetic  guidance  and  contact,  sur¬ 
veillance  of  ground  and  aeaospace  Objects,  Intelligence  data 
collection  and  handling,  Infoamation  system  technology, 
lonosphealc  propagation,  solid  state  sciences,  microwave 
physic*  and  electronic  reliability,  maintainability  and 
compatibility. 


